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ABSTRACT 



We study destriping as a map-making method for temperature- and-polarization data for cosmic microwave background observations. 
We present a particular implementation of destriping and study the residual error in output maps, using simulated data corresponding 
to the 70 GHz channel of the Planck satellite, but assuming idealized detector and beam properties. The relevant residual map 
is the difference between the output map and a binned map obtained from the signal + white noise part of the data stream. For 
destriping it can be divided into six components: unmodeled correlated noise, white noise reference baselines, reference baselines 
of the pixelization noise from the signal, and baseline errors from correlated noise, white noise, and signal. These six components 
contribute differently to the different angular scales in the maps. We derive analytical results for the first three components. This study 
is related to Planck LFI activities. 

Key words, methods: data analysis - cosmology: cosmic microwave background 
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1. Introduction 

Construction of sky maps from the time-ordered data (TOD) is 
an important part of the data analysis of cosmic microwave back- 
ground (CMB) surveys. For large surveys like PlanckQ (Planck 
Collaboration I2005T) . this is a computationally demanding task. 
Methods which aim at finding the optimal minimum- variance 
map (Wright T996', Borrill 1999, Dore et al. 2001, Natoli et 
al. 2001 , Yvon & Mayet 2005, de Gasperis et al. 2005 ) are com- 
putationally heavy and require large computers. Also, a faster 
method is needed for Monte Carlo studies to assess systematic 
effects, noise biases, and error estimates. 

Destriping (Buriganaet al. 1997b, Delabrouille [T998l Maino 
et al. 19M 2002. Revenu et al. 2000 Sbarra et al. [20031 
Poutanen et al. 120041 Keihanen et al. 12004b is a fast map-making 
method that removes correlated low-frequency noise from the 
TOD utilizing crossing points, i.e., the same locations on the 
sky observed at different times. Correlated noise is modeled as a 
sequence of ("uniform") baselines, i.e., constant offsets in the 
TOD. High-frequency noise (frequency of the same order or 
higher than the inverse of the baseline length) cannot be mod- 
eled this way. Thus the method assumes that the high-frequency 
part of the noise is uncorrelated (white noise). 

In some implementations, a set of base functions (e.g. low 
order Legendre polynomials) is used instead of just the uni- 
form baseline (Delabrouille 1998, Maino et al. 2002, Keihanen 
et al. 2004, 2005), or a spline is fitted to the TOD (Ganga 1994). 

In this paper we describe one destriping implementation 
for making temperature and polarization maps of the sky and 
study the residual errors in the maps. This implementation was 



Send offprint requests to: H. Kurki-Suonio, e-mail: 
hannu.kurki-suonio@helsinki . fi 
^ ihttp : //www ■ rssd . esa . int/index . php?pro j ect=PLANCK| 



originally known as the "Polar" code, and used in the map- 
making comparison studies of the Planck CTP Working Group 
(Poutanen et al. l2QQ6l Ashdown et al. 2007 al 12007b] l2QQ9l) . Polar 
has now been merged into the "Madam" destriping code. The 
novel feature in Madam was the introduction of an optional noise 
prior (noise filter) that utilizes prior information on the noise 
power spectrum (Keihanen et al. 2005 ). Polar corresponds to 
Madam with the noise prior turned off. The results presented in 
this paper were obtained with the Madam code, with the noise 
prior turned off. We briefly comment on the effect of the noise 
prior in Sect. [H The use and effect of the noise prior will be 
described in detail in Keihanen et al. (120091) . 

Destriping errors have been previously analyzed by Stompor 
& White (120041) and Efstathiou ( 20051120071) . 

The TOD can be considered as a sum of signal + white noise 
-h correlated ("1//") noise. If there were no correlated noise, 
the optimal way to produce a map would be a simple binning 
of the TOD samples onto map pixels. (We do not address here 
the question of correcting for the effect of the instrument beam. 
"Deconvolution" map-making methods that correct for the effect 
of the beam shape have been developed (Burigana & Saez 120031 
Armitage & Wandelt 2004, Harrison et al. 2008), but tend to 
be computationally very resource intensive. They also alter the 
noise properties of the maps in a way that is difficult to follow 
in CMB angular power spectrum estimation.) Thus the task of 
a map-making method is to remove the correlated noise as well 
as possible, with as little effect on the signal and white noise 
as possible. The difference of the output map from the binned 
signal + white noise map is thus the residual map to consider to 
judge the quality of the output map. We divide this residual into 
six components: unmodeled 1 // noise, 1 // baseline error, white 
noise reference baselines, white noise baseline error, pixelization 
noise reference baselines, and pixelization noise baseline error. 



2 



Kurki-Suonio et al.: Destriping CMB temperature and polarization maps 



We study the nature of each component, and its dependence on 
the basehne length. 

We have used simulated data corresponding to 1 year of ob- 
servations with 4 Planck LFI 70 GHz detectors (two horns, each 
with two orthogonally polarized detectors). 

For simplicity, we did not include foregrounds in the signal 
(see Ashdown et al. 2007b, 2009 for effect of foreground signal) 
or such systematic effects as beam asymmetries, sample integra- 
tion, cooler noise, or pointing errors (see Ashdown et al. 2009 ). 
Even in Ashdown et al. (2009 ) the simulated data used was still 
fairly idealized. We are currently working on more realistic sim- 
ulations. 

In Sect. [2] we discuss some early-stage design choices made 
in the development of our map-making method. Sect. [3] contains 
the derivation and description of the method. Sect. (H describes 
the simulated data used to test the method. In Sect. [5] we an- 
alyze residual errors in the time domain, and in Sect. [6] in the 
map domain. In Sect. [7] we discuss the effect of the noise knee 
frequency, and in Sect. [8] we give a preview of results obtained 
when a noise prior is added to the method. We mainly consider 
maps made from a full year of data, but in Sect. [9] we discuss 
maps made from shorter time segments. In Sect.[TOlwe summa- 
rize our conclusions. 



2. Design choices 

2.1. Ring set or not 

For a PLANCK-like scanning strategy, where the detectors scan 
the same circle on the sky many times before the spin axis of 
the satellite is repointed, an intermediate data structure can be 
introduced between the TOD and the frequency map. The cir- 
cles from one repointing period can be coadded to a ring, i.e., 
averaged to appear just as a single sweep of the circle. In this 
context it is natural to choose one baseline per ring. Destriping 
is then performed on this ring set, instead of the original un- 
coadded TOD. This reduces the memory and computing time 
requirements by a large factor. 

If the scanning is ideal, i.e., the observations (samples) from 
the different circles of the same ring fall on exactly the same 
locations on the sky, destriping coadded rings is equivalent to 
destriping the uncoadded TOD with baseline length equal to the 
repointing period (i.e., one baseline per ring). In this case it is 
also almost equal (for map-making purposes) to destriping the 
uncoadded TOD with baseline length equal to the spin period 
(i.e. one baseline per circle), see Sect. 16.3.21 

In reality, the spin axis will nutate with some small ampli- 
tude, so that the different circles will not scan exactly the same 
path on the sky. The spin rate is also not exactly constant, and 
the detector sampling frequency is not synchronized with the 
spin. Binning the samples first into a ring ("phase binning") 
(van Leeuwen et al. 2002) and then repixelizing the ring pixels 
into map pixels after destriping may then introduce some extra 
smoothing of the data. 

We have chosen to sidestep this intermediate structure and to 
assign the samples directly to map pixels. The baseline length is 
then not necessarily tied to scan circles and rings, and also data 
taken during the repointing maneuvers can be used. For short 
baselines no data compression is possible and map-making is 
done from the full TOD, requiring large computer memory. For 
long baselines (many scan circles) the data can be compressed 
by binning observations directly to map pixels. For baselines one 
repointing period (one ring) long a similar data compression is 
achieved as by using the intermediate ring set structure. Whereas 



the phase-binned ring set is still closely connected to the time 
domain, our "pixel binning" destroys the time-ordered structure 
of the ring, and therefore works only with uniform baselines. 

This way we have achieved a versatile destriping method, 
where the baseline length is an adjustable parameter. Shorter 
baselines can be used in large computers for higher accuracy, 
whereas longer baselines require less memory and computing 
time and can be used in medium- sized computers and for Monte 
Carlo studies. The baseline length is not tied to the scanning 
strategy, and our destriping method can be applied to any scan- 
ning strategy that has crossing points, not just to a PLANCK-like 
scanning strategy. 

However, for a PLANCK-like scanning strategy there is a cer- 
tain advantage in choosing the baseline length so that an integer 
number of baselines fits to one repointing period. Baseline seg- 
ments that extend to two different repointing periods are avoided 
this way. This is mainly an issue for long baselines (not very 
much shorter than the repointing period). For baselines shorter 
than the spin period there seems to be some advantage in choos- 
ing the baseline length so that an integer number of baselines fits 
to one spin period. See Keihanen et al. ( 2009 ). In this paper we 
only consider such choices for baseline length. 

2.2. Crossing points and signal error 

The baselines are estimated from crossing points, i.e., observa- 
tions falling on the same map pixels at different times. Samples 
are assigned to pixels based on the pointing of the detector beam 
center. The beam center may still point at a different location 
within the same map pixel for different samples. We do not at- 
tempt to correct for this effect and this results in a "signal error" 
in our output maps. The signal error due to in-pixel differences in 
beam pointing could be largely eliminated in another kind of de- 
striping implementation, where the scanning circles are treated 
as exact geometrical curves (instead of just a sequence of map 
pixels), and the observations are interpolated to the exact cross- 
ing points of these lines (Revenu et al. 12000b . In this case only 
actual crossings of the scan circles contribute to baseline de- 
termination, whereas in our implementation it is enough that 
two paths pass through the same pixel without actually cross- 
ing there. The latter situation is very common, since successive 
circles are almost parallel. 

However, in a realistic situation there are other contributions 
to signal error that could not be eliminated this way. One such 
contribution is the different beam orientations of the different ob- 
servations of the crossing point, as real beams are not circularly 
symmetric. In Ashdown et al. ( 2009 ) elliptic beams were consid- 
ered for the Planck 30 GHz channel and it was found that this 
had a contribution to the signal error, which was of comparable 
size or larger. 

3. Destriping technique 

3.1. Derivation 

The destriping method can be derived from a maximum- 
likelihood analysis of an idealized model of observations. The 
signal observed by a detector sensitive to one linear polarization 
direction is proportional to 

s = I -\- Qcos2i//-\- Usm2i//, (1) 

where i// is the polarization angle of the detector, and /, Q, and 
U are the Stokes parameters of the radiation coming from the 
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observation direction. In the idealized model the time-ordered 
data (TOD), a vector y consisting of rit samples, is 



y = Pmin + n 



(2) 



where min (the "input map") represents the sky idealized as a 
map of Hp sky pixels, P is the pointing matrix, and n is a vector of 
length Ht representing the detector noise. For observations with 
multiple detectors, the TODs from the individual detectors are 
appended end-to-end to form the full TOD vector y. 

Since we are dealing with polarization data, the map m is 
an object with Srip elements; for each sky pixel the elements 
are the /, Q, and U Stokes parameters. The pointing matrix 
P is of size (rit, 3np). Each row has 3 nonvanishing elements 
(1, cos 20^, sinli//) at the location corresponding to the sky pixel 
in which the detector beam center falls for the sample in ques- 
tion (the sample "hits" the pixel). We do not make any attempt 
at deconvolving the detector beam. Thus the map m represents 
the sky smoothed with the detector beam and the pixel window 
function. The pointing matrix spreads the map into a signal TOD 
Pm. 

We divide the TOD into segments of equal length ^base; 
rif = ^/^^base- For each segment we define an off'set, called base- 
line. The baselines model the low-frequency correlated noise 
component, "1// noise", which we want to remove from the 
data, and we approximate the rest of the noise as white. Thus 
our idealized noise model is 



n = Fain + w 



(3) 



where the vector ain (of length nt) contains the baseline ampli- 
tudes and the matrix F, of size (rit, nt), spreads them into the 
baseline TOD, which contains a (diff'erent) constant value at 
(b = 1, . . . , ^^) for each baseline segment. Each column of F 
contains ^base elements of 1 corresponding to the baseline in 
question, and the rest of the matrix elements are 0. The vector 
w represents white noise, and is assumed to be the result of a 
Gaussian random process, where the diff'erent samples in w are 
uncorrected. 



(4) 



Here ( ) denotes expectation value. The white noise (time- 
domain) covariance matrix 



r 



(ww^) 



(5) 



is thus diagonal (with elements crj), but not necessarily uniform 
(the white noise variance crj may vary from sample to sample). 

Given the TOD y, and assuming we know the white noise 
variance Cw, we want to find the maximum likelihood map niout. 
We assume no prior knowledge of the baseline amplitudes a, 
i.e., they are assigned uniform prior probability. (The variant of 
the method, where such prior knowledge is used, is described in 
Keihanen etal. [2QQ9l ) 

Given the input map m and the baseline amplitudes a, the 
probability of the data y is 



P(y|m,a) = (det27rCJ-i/2 exp(-iw^C;^w) 



(6) 



where w = y - Fa - Pm. This is interpreted as the likelihood 
of a and m, given the data y (Press et al. 1992). Maximizing 
the likelihood is equivalent to minimizing the logarithm of its 
inverse. We obtain the chi- squared function 



X 



-2\nP = (y - Fa - Pm)^ C;^ (y - Fa - Pm) 



(7) 



to be minimized. (We dropped the constant prefactor of Eq. Q.) 
We want to minimize this with respect to both a and m. 

Minimization of Eq. (|7]) with respect to m gives the 
maximum-likelihood map 



m = {v^C-Jvy'v^C-J (y-Fa) 

for a given set of baseline amplitudes a. 

The symmetric non-negative definite matrix 

M = p^c;^p, 



(8) 



(9) 



which operates in the map space, is 3 x 3 block diagonal, one 
block Mp for each pixel p\ 



Mp = 



Z cos 2(/^f 
yi sin l^t 



cos l\f/t 

cos^ 



2. 



sin l\f/t 



2, 

" t 

sin lif/f cos lif/t 



COS 2i//t sin 2i//t 



sin^ l\^t 



(10) 



where the sums run over all samples t that hit pixel p. is the 
white noise covariance matrix for the three Stokes parameters /, 
2, U in pixel p. Mp can only be inverted if the pixel p is sam- 
pled with at least 3 sufficiently diff'erent polarization directions 
i//t, so that all Stokes parameters can be determined. This can be 
gauged by the condition number of M^. If the inverse condition 
number rcond (ratio of smallest to largest eigenvalue) is below 
some predetermined limit, the pixel p is excluded from all maps 
(as are pixels with no hits), and the samples that hit those pix- 
els are ignored in all TODs. Technically, this is done by setting 
M^^ = for such pixels and (C^^)^^ = for the corresponding 
samples. M can then be easily inverted by non-iterative means. 

If all (Tt are equal, crjMp(l, 1), gives the number of hits 
(observations) in pixel p. Thus M is sometimes called the A^obs 
matrix. The optimal distribution of polarization directions i//t 
measured from a pixel is one where they are uniformly dis- 
tributed over 180° (Couchot et al. 1999). In this case = 
(%it,;./cr?)diag(l,l/2,l/2)andM;i = (cr^/^hit,^) diag(l, 2, 2) 
giving the maximum possible value rcond = 0.5. 

Substituting Eq. ([8]) back into Eq. d?]) we get this into the 
form 



;r' = (y-Fa)^Z^C;^Z(y-Fa), 
where we have defined 

z = i-PM-^p^c:/. 



(11) 



(12) 



Here I is the unit matrix. The matrix Z operates in TOD space 
and is a projection matrix, Z^ = Z. If all cr^ are equal, Z is sym- 
metric. In general, C~^Z is symmetric, so that 



c-}z = z' c-/ = z^ C-/Z. 



(13) 



We minimize Eq. (fTTt with respect to a to obtain the 
maximum-likelihood estimate of the baseline amplitudes aout- It 
is the solution of the equation 



(FX'ZF)a = F^C;iZy, 

where we have used Eq. ([T3]) . 
The matrix 

D = F^C:/ZF 



(14) 



(15) 



on the left-hand side of Eq. ([T4l) operates in the baseline space. 
It is symmetric but singular. Eq. ([T4l) has a solution only if its 



Kurki-Suonio et al.: Destriping CMB temperature and polarization maps 



4 

right-hand side is orthogonal to the null space of D. The solution 
becomes unique when we require it to be orthogonal to the null 
space too. 

The null space of D contains the vector that gives all base- 
lines the same amplitude. This represents the inability to detect 
a constant offset of the entire noise stream n, because it has the 
same effect on y as a constant shift in the / of the entire min (the 
monopole). This is of no concern (but should be kept in mind) 
since the goal is to measure the CMB anisotropy and polariza- 
tion, not its mean temperature. If the baselines are sufficiently 
connected by crossing points (two different baseline segments 
of the TOD hitting the same pixel), there are no other kind of 
vectors in the null space, so that the dimension of the null space 
is one. The right-hand side of Eq. (O is orthogonal to this one- 
dimensional null space, and thus Eq. ([T4l) can now be solved. In 
practice it is solved by the conjugate gradient method. If the ini- 
tial guess is orthogonal to the null space, the method converges 
to a solution that is also orthogonal to the null space. Normally 
we start with the zero vector as the initial guess to guarantee this. 
This means that the average of the solved baseline amplitudes is 
zero. Strictly speaking this holds exactly only when no pixels 
are excluded from the baseline determination due to their poor 
rcond. 

We write the solution of Eq. ([T4l) as 

aout = D-iFXiZy. (16) 

D"^ is interpreted as the inverse in this orthogonal subspace. D 
and D"^ will act in this subspace only. 

Using the maximum likelihood baselines from Eq. ([T6l) in 
Eq. ([5]) we get the output map of the destriping method: 

mout = M-ipXHy-Faout). (17) 

Eqs. ([T6l) and ([TTl) summarize the destriping method. 

Implementation details are discussed in Keihanen et 
al. (.2QQ9j . 

3.2. Description 

Let us review the different operations involved: 

P^C^^ acts on a TOD y to produce from it a sum map 
P^C^^y where each pixel has Stokes parameters representing a 
sum over observations that hit the pixel, 

Qp = ^^cos2jA.};, (19) 

tep 

Up = ^^sin2(A,J, (20) 

tep 

Instead of a sum, we should take the average of observa- 
tions. This is accomplished by M"^P^C^^ which corresponds 
to solving the Stokes parameters from the observations hitting 
this pixel; i.e., without regard for pixel-to-pixel noise correla- 
tions. The observations are just weighted by the inverse white 
noise variance. The resulting map is called the naive map, or the 
binned map. We shall use the shorthand notation 

B = M-^P^C;^ . (21) 

B acts on a TOD to produce from it a binned map. Note that 
BP = I. 



(It may be better to use the same at for the two polarization 
directions of the same horn, to avoid polarization artifacts from 
systematic effects (Leahy et al. ' 2QQ9I) . In case their true noise 
levels are different, destriping allows also the option of using 
equal CTt in solving baselines, Eq. ([T6l) , but the actual CTt in the 
final binning to the map, Eq. ([TT]) . We do not study this issue in 
this paper, as we used simulated data with a constant cr^.) 

We can now see that the effect of 

Z = I - PB (22) 

on a TOD is to bin it to a map, read a TOD out of this map, 
and subtract it from the original TOD. Thus Zy represents an 
estimate of the noise part of y. Acting on a TOD constructed 
from a map as Pm, Z returns zero, ZPm = 0. 

Likewise, F^C^^ acts on a TOD to sum up the samples of 
each baseline segment, weighting each sample by cr~^. 

The effect of the matrix 

D = F^C;^ZF = F^C;^F - F^C;^PBF (23) 

on a baseline amplitude vector a is to produce a TOD Fa con- 
taining just these baselines, make a noise estimate ZFa from this 
baseline TOD, and calculate then the weighted sum F^C^^ZF of 
this noise estimate for each baseline. Thus the content of Eq. ([T4l) 
is to find such a set of baseline amplitudes that these noise es- 
timate sums are the same for the baseline TOD Fa as for the 
actual input TOD y. Thus the solved baselines aout represent the 
estimated average noise for each baseline segment of the TOD. 

These baselines are then subtracted from the TOD to produce 
the cleaned TOD y - Faout, which is then binned to produce the 
output map 

mout = B(y-Faout) (24) 
in Eq. (fTTl) . 

In a good scanning of the sky the number of hits in each pixel 
is large. From Eq. ([22b we see that Z contains two parts. The first 
part I gives each row t a large diagonal element 1 corresponding 
to the TOD sample this row is acting on. The second part gives 
this row a large number of small nonzero elements correspond- 
ing to all samples f that hit this same pixel. The sum of these 
elements is -1 so that the row sum is zero. Thus the first and 
second parts make an equally large contribution, but the second 
part comes in many small pieces. 

The matrix D has a similar structure. The first part (see 
Eq. [23]) is diagonal containing the sum Yuteb 1 l^li ^^^^ 
pies in the baseline segment b. The second part gives to each 
row b a nonzero element for each baseline b^ that has a cross- 
ing point with b. For a good scanning each baseline has a large 
number of crossing points, so that this second part contributes a 
large number of small elements to each row. 

We define a shorthand notation for the matrix 

A = D-^F^C;^Z, (25) 

which appears in Eq. ([T6b as aout = Ay. This matrix acts on 
a TOD and produces from it the set of baseline amplitudes ac- 
cording to the destriping method, Eq. ([T6b . Note that AF = I, in 
the sense that AFa = a for any set of baseline amplitudes, whose 
average is zero (so that a is orthogonal to the null space). Since 
Z is a projection matrix, AZ = A. 

Note that all the matrices introduced (see Table [B are con- 
structed from I, P, Cw, and F. The adjustable parameter in the 
destriping method is the baseline length ^base, which affects the 
matrix F. P is determined by the scanning strategy and map pix- 
elization, and Cw by detector noise properties. 
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01 
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pointing matrix 


F 


Uf X 




UCl5>dillC5> IVJ 1 KJLJ 


Cvv = (ww^> 


Ht X 


121) 


white noise cov. 


M = P^C; P 


3;ip X 3np 


19]) 


A^obs matrix 


B = M ip^C;i 


3np X 


IID 


bin TOD to a map 


Z = I - PB 


X Ht 


IE]) 




D = F^C; ZF 


ntxrib 


123 




A = D F^C; Z 


fib X fit 


125]) 


solve baselines 


R = (f^C;1f)"' F^C;1 


fib X 


07]) 


reference baselines 



Table 1. Table of matrices. The square matrices C^, M, D, and 
C^^Z are symmetric, is diagonal, and M 3x3 block diagonal. 
Z is a projection matrix. D is singular, and D"^ is its inverse in 
the subspace orthogonal to its null space. The third column refers 
to the equation in which the matrix was introduced. 
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3.3. Destriping error 
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i.e, the matrix D of Eq. ([T5]) is the Fisher matrix of the baselines. 
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Its inverse gives the covariance of the baseline error Aa = 
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Ignoring the nondiagonal terms we get an approximation 



Multipole, f 
EE 



(28) 



Assuming the white noise variance stays constant, at = cr, over 
a given baseline, this becomes 



'^base 



(29) 




We can understand the approximate result (l29l) as follows. 
Destriping solves the baseline amplitude at from the differences 
between samples from baseline segment b and from other base- 
line segments that hit the same pixel. Both the baselines ain 
and the white noise contribute to these differences. Without the 
presence of white noise the baselines could be solved exactly, 
^out = ^in (in the idealized model considered in this section), so 
the error comes from the white noise in the samples. Typically 
there are many more differences involved than there are sam- 
ples in the baseline (number of hits per pixel » 2), but samples 
from the baseline segment b contribute only ^base uncorrected 
random variables, with variance cr^, and their contribution cor- 
responds to the result (|29b . The total of samples from baseline 
segment b contribute equally to its amplitude determination than 
do the total of samples from other baseline segments that hit 
the same pixels. But the latter contribution comes from a much 
larger number of uncorrected random variables, with a similar 
variance, so it gets averaged out by the larger number and makes 
a much smaller contribution to the baseline error variance dZTl) . 
This contribution from the white noise from the crossing sam- 
ples is ignored in (|28l) and (|29l ). However, it should add to the 
baseline error variance, so that Eqs. (l28l) and ([29b are underesti- 
mates. See also Efstathiou (l2QQ5ll2QQ7l) . 



Fig. 1. The input angular power spectrum. Black lines show the 
theoretical Q produced by CAMB, and the beam-and-pixel- 
smoothed version of it. Red dots show the Q of the a^m real- 
ization. Blue dots show the Q calculated from the binned noise- 
less map. (A version of this paper with higher-quality figures is 
available at http://www.helsinki.fi/~tfo_cosm/tfo-planck.html ) 



The approximation (|29l ) corresponds to the white noise refer- 
ence baseline contribution discussed in Sect. [5] We see in Sect. [6] 
that, while this approximation is good in the time domain, it is 
not that good in the map domain. 



4. Simulation 

In this study we tested the destriping method using simulated 
data that is more realistic than the model used to derive the 
method in Sect. [3] The TOD was prod uced u sing Planck Level- 
S simulation software (Reinecke et al. I20Q6I) as a sum of signal, 
white noise, and correlated noise (called also 1 // noise). 



(30) 
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Fig. 2. Binned CMB signal map, / and Q, full sky. 
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Fig. 3. Binned CMB temperature (/) signal maps from two 10° x 
10° regions. Left column: Near the ecliptic North Pole (centered 
at ^ = 7°, = -90°). Right column: Near the ecliptic (centered 
at6/= 85°,0 = 5°). 



We produced the three time streams, s, w, and separately, so 
we could analyze the effect of each component on the destriping 
error. The TOD represented 1 full year (366 d) of data from 4 
polarized Planck 70 GHz detectors (LFI-19a, LFI-19b, LFI-22a, 
LFI-22b). 



4.1. Signal 

We considered the CMB signal only; no foregrounds were in- 
cluded in the simulation. Detector pointings 6t,(pt,^t for each 
sample t were produced to imitate a realistic scanning strategy. 
We used a set of input spherical harmonic coefficients a^^, a^^ 
to represent the sky. The signal sample St was then produced by 
the convolution of a circularly symmetric Gaussian beam (fwhm 
= 12.680 centered at this pointing with the input a^^ (Wandelt 
& G6rski,200L Challinor2Q0a). Thus the /, 2, and U (Eq.\B 
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Fig. 4. Same as Fig. [3l but for the Stokes parameters Q (top) 
and U (bottom). Since the signal contains only E mode polar- 
ization, Q shows structures elongated along lines of latitude and 
longitude, whereas U shows structures elongated 45° away from 
them. 



of the signal parts of different samples hitting the same pixel are 
different, as Ot, (pt can vary within the pixel. 

To produce the input we used the CAMBEI code to pro- 
duce the theoretical angular power spectra Cj^, Cj^, C^^ for the 
Friedmann-Robertson- Walker universe with cosmological pa- 
rameter values Qo = 1, = 0.7, cOm = 0.147, cob = 0.022, 
T = 0.1, and with scale-invariant (^ = 1) adiabatic primordial 
scalar perturbations with amplitude 5 x 10"^ for the curvature 
perturbations. A realization was then produced from 

these spectra. Effects of gravitational lensing were ignored, and 
therefore there is no B mode polarization in the input spectrum. 

Fig. [T] shows the input Q as well as the Q of the binned 
(noiseless) signal map Bs made from the simulated TOD. 

Fig. [2] shows the binned signal maps Bs, and Figs. [3] and (H 
zoom into two 10° x 10° regions to reveal small-scale detail. In 
this paper we keep using these same two 10° x 10° regions, one 
near the ecliptic north pole, one near the ecliptic, to show map 
detail. 



4.2. Scanning 

The sky scanning strategy was cycloidal (Dupac & Tauber 
120051) : the satellite spin axis was repointed at 1 h intervals, caus- 
ing it to make a clockwise circle (radius 7.5°) around the anti- 
Sun direction in 6 months. 

Random errors (rms = 1.3') were added to the repointing. 
Between repointings the spin axis nutated at an amplitude re- 
lated to the repointing error. The mean nutation amplitude was 
1.6'. The nutation was dominated by a combination of two pe- 
riods, ~ 45 s and ~ 90 s. In reality, the repointing errors and 



http://camb.info 
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Fig. 5. Hit map for the 1-year simulation. We show regions 
around the ecliptic North Pole (left) and South Pole (right). The 
color scale is linear and goes from zero (blue) to 50 000 (red). 
Lines of latitude are drawn at 5° intervals. 
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Fig. 7. Hit maps of the two 10° x 10° regions. The color scale 
is linear and goes from to 50 000 in the left plot, and from 
1000 to 3000 in the right plot. The 1-year survey begins and 
ends near the right edge of the right plot, where the first and last 
rings overlap to produce a higher hit count. 
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Fig. 6. Power spectrum of the noise stream n (black). The spec- 
trum was calculated from the detector 19a noise stream by divid- 
ing it into 6-day pieces, obtaining their spectra separately, tak- 
ing their average, and binning the spectrum into 800 logarithmic 
bins, to get a smooth curve. The red curve is the noise model. 



at the four corners of these caustics. For destriping, such a clus- 
tering, where very many crossing points fall on the same pixel, 
is disadvantageous, since there is less independent information 
available for solving the baselines of these rings. The clustering 
occurs when the curvature of the path of satellite pointing on 
the sky equals the curvature of the scanning circle. The curva- 
ture changes sign when the spin axis is close to the ecliptic, but 
slightly north of it, and the clusterings near the north and south 
ecliptic poles occur a little bit before and after that. Conversely, 
the crossing points are spread more widely along the caustics 
when the spin axis is near its north or south extrema. 

The sampling frequency was /sample = 76.8 Hz, and each 
sample was simulated as an instantaneous measurement (no in- 
tegration along scan direction). The baseline length ^base, is 
given in time units as ^base = ^base //sample in the following. The 
beam center moves on the sky at an angular velocity ^^^scan = 
360°/spin sin ^det ~ 6°/s, so that one baseline corresponds to a 
path of length ^base = ^scan^base on the sky. The samples are sep- 
arated by 6s = tLJscan //sample ~ 4.68' on the sky. The length of 
the simulated TOD was rit = 4 x 366 x 24 x 3600 s x /ample = 
9714401280. 



nutation amplitudes are expected to be smaller. Thus eff'ects of 
nutation and small-scale variations in the map pixel hit count 
appear somewhat exaggerated in this study. 

The satellite rotated clockwise (i.e., spin vector pointing 
away from the Sun) at about 1 rpm (fpm ~ 1/60 Hz); spin 
rate variations (rms 0. 1 °/s) around this nominal rate were chosen 
randomly at each repointing. (In reality, the spin rate variations 
are expected to be smaller.) Coupled with the 60 s spin period, 
the 45 s and 90 s nutation periods produce a 3-min periodicity in 
the detector scanning pattern. 

We simulated 4 detectors corresponding to 2 horns (19 and 
22). The detectors were pointed ^det = 87.77° away from the 
spin axis, causing them to draw almost great circles on the sky, 
the "22" trailing the "19" detectors by 3.1° , but following the 
same path. The a and b detectors of each horn shared the same 
pointing but had diff'erent polarization directions by exactly 90°. 
The polarization directions of the 19 and 22 detector pairs dif- 
fered from each other by 44.8°. 

The cycloidal scanning causes the detector scanning rings to 
form caustics around the ecliptic poles, where nearby scanning 
rings cross (see Fig. [5]). A large number of ring crossings cluster 



4.3. Noise 

The noise part was produced as a sum of white and correlated 
noise, 

n = w-hnc (31) 

where the correlated part (1// noise) was produced by a 
stochastic-diff'erential-equation (SDE) method, that produces 
noise whose power spectrum is approximately of the form 
Pc(f) oc f~^. It is not of the form Fa, but contains a part that 
cannot be modeled with baselines. 

The white noise rms was set to cr = 2700 yuK (thermody- 
namic scale for CMB anisotropics), corresponding to the Planck 
70 GHz goal sensitivity (Planck Collaboration 120050 . The 1// 
noise was simulated with slope -Of = -1.7and/niin = 1.15x10"^ 
Hz (period of one day), so that the power spectrum was flat for 
/ < /min- See Fig.[6j Since the white and 1// streams were pro- 
duced separately, the knee frequency fk (where the white and 
1 // noise powers are equal) could be adjusted by multiplying 
the 1 // stream with difl'erent factors. We used fk = 50 mHz as 
the reference case, representing a conservative upper limit for 
the Planck 70 GHz detectors (Burigana et al. 1997ai Seifl'ert et 
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Fig. 8. Destriped (output) temperature map (one year survey, 15 
s baselines) for the two 10° x 10° regions. 



al. 120021 Tuovinen 120031) . but consider also fk = 25 mHz in 
Sects.[7]and[8]. 

The power spectrum of the 1 // noise is thus approximately 



Pcif) = 




where fc = /sample /2 = 38.4 Hz is the Nyquist (critical) fre- 
quency. 

The l/f noise was generated in ^gen = 15.25 d pieces. 
The actual statistics calculated from the simulated 1 // stream 
were: mean = -9.86yuK, stdev = (Tc = 1868.5yuK, so that 
cr^ = 0.48(T^. As can be seen from Fig. [6l the simulated l/f 
stream power falls below the noise model for the lowest frequen- 
cies. A sizable part of the variance in the 1 // stream comes from 
the very lowest frequencies. 

The nonzero mean of the 1 // stream has to be taken into ac- 
count when comparing solved baselines to the input 1 // stream, 
since the destriping method sets the average of the solved base- 
lines to zero. 

4.4. Maps 

The maps were produced in the HEALPi^Jl pixelization (Gorski 
et al. 2005) in ecliptic coordinates. We used the A^side = 512 
resolution for all maps, corresponding to = 3 145 728 for the 

1 /2 

full sky, with square root of pixel solid angle Op = flp = 6.87'. 
For the full-year TOD the polarization of each pixel was well 
sampled; the lowest rcond was 0.422. The mean rcond was 0.492 
and the maximum 0.49999. The hit count ^hit (number of hits per 
pixel) varied from 818 to 273480 for the full 1-year simulation. 
The mean number of hits was (%it) = rit/rip = 3088.125. The 
mean inverse hit count was (^j^-J) = 0.000418 = (0.0204)^ ^ 
1/2391. Fig. \5\ shows the hit count in the regions around the 
ecliptic poles, where it varies a lot. Fig. [7] shows the hit count in 
the two 10° X 10° regions. 

Fig. [8] shows the output temperature map niout for the case 
^base = 15 s. The visual appearance is the same for other baseline 
lengths. To see differences one has to look at residual maps, see 
Sect.[6l Polarization maps (Fig.O are dominated by small-scale 
noise. The most obvious map-making related feature in the out- 
put maps is the reduction of noise where the hit count is larger. 
Other effects are more subtle and are analyzed in the following 
sections. 




Fig. 9. Even near the ecliptic poles, where the noise in the output 
map is the lowest, the pixel-scale noise from four 70 GHz detec- 
tors is higher than the CMB polarization signal in the map. Left: 
Binned signal Q map. Right: Output Q map (15 s baselines). 
(Both plots are from the same pole region.) 



5. Time domain 

We now analyze the application of the destriping method to this 
kind of data. We consider the case of the full 1-year survey with 
fk = 50 mHz noise, except for Sect. [71 where we discuss the 
effect of changing the knee frequency, and for Sect. [9l where we 
consider maps made from shorter pieces of the TOD. Consider 
this first in the time domain, i.e., look at the cleaned TOD 

d = y-Faout. (33) 

We assume that the two noise streams, w and are indepen- 
dent Gaussian random processes. 

Since the destriping method is linear, we can divide the base- 
line amplitudes obtained by Eq. ([T6l) into the parts coming from 
the different TOD components, 

Rout = Ay = As -h Aw -\- Auc . (34) 

Likewise, the cleaned TOD can be divided into five terms, 

d = s -h w - FAw -h (Uc - FAiic) - FAs . (35) 

The first term is the signal and the second term is the white noise. 
We call the third term white noise baselines, the fourth term (in 
parenthesis) residual 1 // noise, and the fifth term signal base- 
lines. 

The TOD vector Zs appearing in the signal baseline term 

FAs = FD-^F^C;^Zs (36) 

is the pixelization noise (Dore et al. 1200 lb . It is the noise estimate 
we get from the signal TOD (which contains no noise). Signal 
gradients (of the beam- smoothed input sky) within a map pixel 
are the origin of the pixelization noise. 

5.1. Reference baselines 

We define the reference baselines Rn of a noise stream n as the 
weighted averages of each baseline segment, i.e., matrix R is 
defined as 

R^(fX'f)"'f^C;^ (37) 

Note that RF = I. 

The reference baselines of the 1 // noise. 



^ http://healpix.jpl.nasa.gov 



2ir = Rile , 



(38) 
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can be viewed as the "goal" of baseline estimation. Subtracting 
them from the full TOD gives us a TOD stream 

y-Fa, = s + w + (n,-Fa,) , (39) 

which contains, beside the signal and the white noise, only the 
part (lie - Fa^) of correlated noise that cannot be represented in 
terms of baselines. We call this unmodeled 11 f noise. 

The actual cleaned TOD that results from destriping, can 
now be written as 

d = y-Faout = s+w-FAw+(nc - Fa,)-F (An^ - a,)-FAs , (40) 

where the residual 1// noise is split into the unmodeled 1// 
noise and the 1// baseline error, An^ - a^ = (A - R)nc. Thus, 
altogether, there are three contributions to missing the goal of 
baseline estimation, 

^out - = Aw + (A - R) He + As , (41) 

white noise baselines, 1 // baseline error, and signal baselines. 

Unlike the solved baseline contributions An, the reference 
baselines Rn do not involve the pointing matrix (except for the 
case of pixelization noise), so the differences between them are 
related to how the scanning strategy connects the baselines with 
crossing points. Thus, for analyzing errors, it is useful to sep- 
arate also the white noise baselines into white noise reference 
baselines and white noise baseline error. Aw = Rw + (A - R)w; 
and likewise the signal baselines into reference baselines of pix- 
elization noise and signal baseline error As = RZs -h (A - R)Zs. 

The white noise baseline error stream (A - R)w is uncor- 
rected with the white noise reference baseline stream Rw. To 
show this, we note that 

RCwA^ = RC^(l -B^P^)c; FD^ 
= D"^ - RPBFD"^ 

= - (f^Cf)"' (f^Cf - d)d-i 

= (f^Cf)"' = RC,R^ (42) 

so that 

((Rw) ((A - R)w)^) = RC^A^ - RCwR^ = . (43) 

5.2. Approximation to soived baseiines 

The solved baselines aout can now be written as 

aout = ar + As -h Aw -h A(nc - Far) 

+ (f^c;^zf)"'f^c;1zw + 

+ (f^C;^ZF)"' F^C;^Z(n, - Fa,) (44) 
(up to an overall constant), so that 

F^C;^ZF(aout - a,) = F^C;^Zs + F^C^^Zw + 

+ FXiZ(n,-Fa,). (45) 

For the TOD streams F(aout - a,), w, and - Fa, there should 
be no significant correlations between the samples from diff'erent 
circles that hit the same pixel. (If the baseline length is longer 
than one scanning circle, this statement is limited to the samples 
that come from diff'erent baseline segments for F(aout - ^r) and 
lie - Fa,.) Therefore, for a large number of hits, the eff'ect of the 
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PB part of Z tends to average out, leaving the I part dominant, 
and we can approximate Z ^ I, leading to A R, for these 
terms. We get the approximation 

aout ~ a, + Rw -h RZs . (46) 

Comparing to Eq. (|4T]) . we note that in Eq. (l46b white noise base- 
lines are approximated by white noise reference baselines, 1 // 
baseline error is ignored, and signal baselines are approximated 
by the reference baselines of pixelization noise. Thus, for a good 
scanning, the solved baselines aout should track the reference 
baselines of instrument noise -i- pixelization noise. 

5.3. Division 

The eff'ect of destriping on the white noise is just harmful for 
the maps (we elaborate on this in Sect. O, so the relevant time 
domain residual is 

d-s-w = lie - Faout 

= -FAw -h (lie - FAiie) - FAs . (47) 

It consists of three components: 

1. white noise baselines, 

2. residual 1 // noise, 

3. and signal baselines, 

which are uncorrelated with each other. 

Each component can be further divided into two parts: 

d-s-w = -FRw-F(A-R)w 

-h(iic - Fa^) - F(Aiic - a^) 

-FRZs - F(A - RZ)s . (48) 

We call these six components: 

la) white noise reference baselines 

IjS) white noise baseline error 

2a) unmodeled 1 // noise 

213) 1 // baseline error 

3a) reference baselines of pixelization noise 

313) signal baseline error 

Of these components, la, 1(3, 2 and 3 are uncorrelated with 
each other. This division forms the basis of the discussion in the 
rest of this paper. Approximation (l46l) corresponds to ignoring 
the yS-components. Components la and 2a are independent of 
scanning strategy, which makes their properties easier to under- 
stand. The yS-components are related to how the baselines are 
solved using crossing points and thus couple to the scanning 
strategy in a complicated manner. 

We turn now to our results with simulated data to see how 
this comes out in practice, first for the noise part, and then for 
the signal part. 

In Fig. [TOl we show a part of the simulated (1// + white) 
noise stream n, both the reference and solved baselines, FRn and 
FAn, and their diff'erence F(R - A)ii = F(A - R)iic + F(A - R)w, 

for ^base = 15 s. 

We consider now separately the white noise and 1 // noise 
parts. 

5.4. White noise baseiines 

From Tables [2] and [3] we see that the white noise baselines track 
their reference baselines well. The white noise reference base- 
lines are just white noise themselves, their variance cr^^ down 
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Fig. 10. First 15 minutes of the (1// + white) noise stream 
(grey) and its reference (black solid) and solved (black dashed) 
^base = 15 s baselines. Note that these two baseline curves are 
nearly on top of each other. The difference of the curves is plot- 
ted as the red curve and corresponds to the scale on the right. 
This represents the error in the approximation (|46l) , except for 
the signal contribution to it. 
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Fig. 11. Autocorrelation function for the white noise baseline 
error (A - R)w for ^base = 1 hour (black), 1 min (red), 15 s 
(green), 2.5 s (blue). The solid lines are for the "temperature" 
baselines (a-\-b)/2 and the dashed lines are for the "polarization" 
baselines (a - b)/2. For the ^base = 15 s and 2.5 s cases, only lags 
that are multiples of 1 min (^ the spin period) are included. 



^base 


Eq. 


(Twr (reference) 


cr^b (solved) 


2.5 s 


194.856 


194.845 


195.115 


15 s 


79.550 


79.551 


79.631 


1 min 


39.775 


39.766 


39.808 


Ih 


5.135 


5.130 


5.399 



Table 2. Statistics of the white noise baselines. We give the 
standard deviation (square root of the variance) of the baselines 
(in yuK). 



Fig. 12. Autocorrelation function for the white noise baseline 
error (A - R)w for ^base = 2.5 s and detector 19a. For short 
lags, only baselines whose lag is a multiple of 1 min (^ the spin 
period) show significant correlation with each other. For longer 
lags this distinction disappears due to random spin rate varia- 
tions. Since these variations have a rms which is about 1 / 60 of 
the spin rate, for lags around 30 min any one of the baseline seg- 
ments within a spin period is about equally likely to land on a 
given location of the scan circle. The correlation between two 
baselines that land on the same location of the circle is much 
larger (presumably similar to the ^base = 1 min and 1 hour cases), 
but the way we calculate the autocorrelation function (in time 
domain, not in the spin phase domain) is not able to pick this 
out. 



^base 


a and b 


(a + b)/2 


(a - b 


2.5 s 


10.244 


5.916 


8.354 


15 s 


3.595 


2.063 


2.944 


1 min 


1.854 


1.066 


1.517 


Ih 


1.719 


0.989 


1.406 



Table 3. Statistics of the white noise baseline error. We give the 
standard deviation 0-^^ in yuK. 



^base 


(a + b)/2 


(a - b)/2 


overlap 


2.5 s 


0.826 


-0.003 


0.793 


15 s 


0.901 


0.005 


0.966 


1 min 


0.925 


-0.011 


0.991 


Ih 


0.998 


-0.013 


0.99986 



Table 4. Correlations between detectors 19 and 22 for white 
noise baseline error. 



from the white noise variance cr^ by ^base- The white noise base- 
line variance 

-^.--^.~ = (f)-^ = 7^> (49) 

'^tbase \Jc / J5%ase 

is slightly larger. 

In Table [3] we show the statistics of white noise baseline er- 
ror. We show it also for the average and the diff'erence between 
the two polarization directions a and b, which represent the con- 
tribution of this eff'ect to the temperature and polarization mea- 
surements. 
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The white noise reference baselines are completely uncor- 
related with each other. This is not true for the solved white 
noise baselines. The difference, the baseline errors, show sig- 
nificant autocorrelation, see Figs. [TT] and [121 and correlation be- 
tween detectors 19 and 22, see TablelH Although the white noise 
baseline variance is not much larger than the white noise refer- 
ence baseline variance, these correlations make the diff'erence 
between them important. 

These correlation properties are easy to understand. While 
the amplitude of the reference baseline arises from the noise of 
the baseline segment itself, the error (A - R)w is caused by the 
noise in the crossing baselines. Since the baseline segments that 
are separated from each other by an integer number of spin pe- 
riods, and not too many pointing periods, cross almost the same 
set of other baseline segments, often in the same pixels, their 
baseline errors (A - R)w are strongly correlated with each other. 
Also, since the corresponding baseline segments from the horns 
19 and 22 are only 6i;scan3.1° = 0.517 s shifted from each other, 
they crossed almost the same set of other baseline segments, and 
are thus strongly correlated. The overlap fraction 1 -0.517 s/^base 
is given in Table (H Note that this number does not take into ac- 
count that the hits from the two horns may still be distributed 
diff'erently to the pixels in the overlap region. 

More exactly, the "temperature" combinations (a -h b)/2 of 
(A-R)w are strongly correlated between 19 and 22, whereas the 
"polarization" combinations (a - b)/2 are not. The polarization 
directions of the ab pairs 19 and 22 diff'er from each other by 
44.8° making their polarization measurements (a - b)/2 almost 
orthogonal. Thus they also pick almost orthogonal error com- 
binations from the crossing baseline segments, and remain un- 
corrected. This also explains why the stdev of the "polarization 
baseline error" is larger than the "temperature" one in Table [3] 
In eff'ect, only half of the crossing baseline pairs ab contribute to 
determining an (a - b)/2 baseline combination, so the number of 
degrees of freedom is down by 1 / 2 and the variance thus larger 
by 2. 

5.5. Residual 1/f noise. 

Fig. [13] shows the l/f part of the noise n^, its reference base- 
lines sir and the solved l/f baselines, An^. The diff'erence be- 
tween these two sets of baselines, 1 // baseline error, is shown 
in Fig.[T4l 

The 1 // noise can be separated into reference baselines and 
unmodeled l/f noise, = Fa^ + (n^ - Fa^). When we sep- 
arate the 1 // baseline error, (A - R)nc into the corresponding 
components, 

(A - R)nc = (A - R)Fa, -h A(nc - Fa,) , (50) 

we note that the first term on the right hand side gives the same 
contribution to each baseline (= - the mean of the l/f noise), 
and is thus irrelevant. Thus the 1 // baseline error arises from 
the unmodeled l/f noise, in the same manner as white noise 
baseline error (A - R)w arises from white noise. (The reference 
baselines of unmodeled 1 // noise are zero). 



5.5.1. Unmodeled 1// noise 

Since baselines can only model frequencies < fx = l/(2^base) 
and do it better for lower frequencies, the power spectrum of the 
unmodeled 1 // noise - Fa, is equal to that of for / » /x 
and falls rapidly towards smaller / for / < f^. See Fig.[T5l 

Since the reference baselines a, are obtained from the stream 
He, the power spectrum Puif) of the unmodeled l/f noise 




1 /f noise 

Solved baselines 

Reference baselines 



200 400 
time [s] 



Fig. 13. First 10 minutes of the l/f noise stream together 
with its solved and reference baselines for ^base = 1 min. The 
diff'erence between the 1 // noise and its reference baselines is 
the unmodeled 1 // noise. 

He - Far can be obtained from the power spectrum Pdf) of 
through a transfer function. 



Pu(f) = H(f)P,(f) . 
This transfer function is 



H(f) 



I 



sin^(7r/rbase)'^ 



(TrAase)^ 



(51) 



(52) 



For / « fx, H(f) oc /4, so the slope of Pu(f) is 4 - for / « /;,. 

We get a rough estimate of the total power in - Fa^ by in- 
tegrating the original Pdf) from to fc and (f/fxfPdf) from 

..^nin to fx. 



I 



5 - a 



a-l 
I 



1-1^ 



a-l 



For a = I this becomes 

ifk 



fc 



2fk(fk 
fc \f. 



a-l 



I 



a-l 



^ -1 



4 fk 



(53) 



(54) 



For our case, a= \.l,a = 2700juK, and we can approximate 

\a-l 

^ 1 , (55) 



-(f 



SO that Eq. ([53l) gives 
cr, ^57.3yuKw(s)'-'^ 



oc t 



ia-l)l2 
base 



(56) 

In Table [5] we compare this approximation to a numerical inte- 
gration of Eq. (|5T1) and the actual variance of - Fa^ from our 
simulated data. We see that Eq. ([56b is an underestimate. An ex- 
ception to this is the 1 -hour-baseline case; this is explained by 
the simulated noise having less power at low frequencies than 
the noise model. 
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Fig. 14. The 1 // baseline error, i.e., the difference between the 
solved and reference 1 // baselines FAn^ - Fa^. The top panel 
shows the case for ^base = Ih (red), 15 min (green), and 4 min 
(blue) from the first 10 hours and from another set of 10 hours 
from the second month of the simulation. The bottom panel 
shows the case for ^base = Imin (red), 15 s (green), and 5 s (blue) 
from the first 10 minutes and from another set of 10 minutes 
from the second month of the simulation. The horizontal line 
at 9.86yuK shows the mean diff'erence that is due to the nonzero 
mean of the 1 // noise. Only the deviations from this mean are 
significant. 



^base 


Eq. 1561) 


Eq. ISB 


actual 


7.5 s 


115.9 


103.1 


121.7 


15 s 


147.8 


131.7 


154.1 


1 min 


240.0 


214.4 


250.3 


Ih 


1006.0 


897.3 


938.0 



Table 5. The standard deviation cr^ (square root of the variance) 
of the unmodeled 1 // stream. Units are yuK. The third column 
is from a numerical integration of Eq. (|5T1) . The last column is 
the actual cTu calculated from the first 16 hours of the simulated 
LFI-19al// stream. 



LFI-19a, 1 min 




Frequency [Hz] 

Fig. 15. Power spectrum Pc(f) of the 1 // noise (black) and its 
unmodeled part Pu(f) (red), for ^base = 1 min. The blue line is an 
analytical estimate corresponding to Eqs. (1321) . ([5T1) . and (1521) . 

1/f noise baseline nnlsestimatlon 




0.0 



-0.2 I \ , , , \ , , , \ , , , \ , , , \ , , , 

2 4 6 8 10 

time [h] 

Fig. 16. Same as Fig. [U but now for the 1 // noise baseline 
error. 



5.5.2. 1// baseline error 

Since the unmodeled 1 // noise is correlated, the properties of 
1 // baseline errors diff'er from white noise baseline errors. Over 
timescales <$c ^base the unmodeled 1 // is positively correlated, 
but since it averages to zero over each baseline segment, there is 
a net anticorrelation. 

From Fig. [14] we see for short baselines a spin- synchronous 
pattern. For baselines of ^base = 1 min, the remaining pattern 
shows a 3-minute period that comes from spin-axis nutation. For 
even longer baselines the long-time- scale correlation shows up 
clearly, l/f baseline error decreases for longer baselines, but 
much less steeply than white noise baseline error. See Table (6] 
Fig. [16] shows the l/f baseline error autocorrelation. For white 
noise baseline error we noted that we get a correlation for nearby 
baselines since they often cross the same other baseline in the 
same pixel. Since unmodeled 1 // noise is positively correlated 
over short timescales, it is not necessary for this crossing to oc- 
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Fig. 17. The 1 // baseline error, FAn^ -Fa^, for 1-hour baselines 
for the full 1-year simulation. We show the sum (red) and the 
difference (blue) between the two polarization directions. The 
times when ring crossings occur near the corners of the caustics 
around the ecliptic poles are centered roughly around hours 240, 
1950, 4610, and 6340. Correlations in baseline amplitudes are 
enhanced near these times. 
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Fig. 18. Effect of baseline subtraction on the noise power spec- 
trum. The solid black line is the spectrum of the original noise 
stream n. The colored lines show the spectra of cleaned noise 
streams for different baseline lengths. Solid lines are for the case 
of subtracting the solved baselines, n-FAn, and the dashed lines 
for subtracting reference baselines, n - FRn. 



^base 


a and h 


(a + h)ll 


(a - b)l2 


5s 


1.790 


0.995 


1.486 


15 s 


1.894 


1.075 


1.557 


1 min 


1.404 


0.765 


1.177 


1 h 


1.360 


0.726 


1.142 



Table 6. Statistics of the 1 // baseline error. We give the stan- 
dard deviation (Tub of the baseline difference An^ - Rn^ in jjK. 



cur in the same pixel to get the positive correlation. Therefore 
we now get positive correlations for even longer timescales. This 
makes 1 // baseline error more important than its small variance 
suggests, since it is not averaged away when binned onto the 
map. 

In Fig. [T7] we show the 1 // baseline error for ^base = 1 h 
over the full year of the simulation for the detector pair 19. The 
effect of the 6-month period of the cycloidal scanning is clearly 
visible. 

5.6. Noise power spectra 

We show the power spectra of the cleaned (destriped) noise 
streams n-FAn for different baseline lengths in Fig.[T8l They are 
compared to the spectrum of the original noise stream n and the 
spectra where noise reference baselines are subtracted instead, 
n - FRn. 

Subtraction of baselines suppresses noise at / < 1/^base- For 
lower frequencies the noise is suppressed more, as baselines can 
model the lower frequencies better. When reference baselines are 
subtracted, the noise power keeps going down toward lower fre- 
quencies; however the solved baselines seem to be able to sup- 
press noise power about 6 orders of magnitude only. 

For shorter baselines, spectral features appear at special fre- 
quencies. They do not appear when reference baselines are sub- 
tracted, so they are clearly related to the scanning strategy. 1) 
There are peaks at / = 1/min and / = 1/(3 min) corresponding to 



the spin and nutation frequencies, and their harmonics. 2) There 
is a notch in power at / = 1/h, corresponding to the repointing 
period, and its harmonics. These are easy to understand: 

The solved baselines come from a noise estimate based on 
subtracting from each sample the average of all samples that hit 
the same pixel. Consider an ideal scanning where the same pixel 
sequence is hit during each spin period within a repointing pe- 
riod: 

1) If the spin period is equal to or a multiple of the period of 
a noise frequency component, all samples hitting the same pixel 
during a given repointing period get the same value from this 
noise component. Thus this noise component can be detected 
as noise (and not signal) only by comparing hits from different 
repointing periods, resulting in a much poorer noise estimate. 

2) On the other hand, if the repointing period is equal to or 
a multiple of the noise period, but the spin period is not, then 
the different samples hitting the same pixel average to zero. This 
noise component is then recognized as noise in its entirety, and 
the solved baseline becomes equal to the reference baseline. 

Since in our simulation the scanning deviated from ideal, 
these features are weakened, but still clearly visible. 

As mentioned in Sect. [T] and Sect. 15.31 and elaborated in 
Sect. [6j the relevant residual noise is the stream n^ - FAn where 
both white noise and l/f baselines are subtracted from the 
1 // noise stream. This is shown in Fig. \19\ The subtraction of 
the white noise baselines has added power to the cleaned 1 // 
stream. At low frequencies the spectrum appears now white, 
since the white noise reference baseline stream is white at 
timescales longer than ^base- The power of the baseline error rises 
towards the lowest frequencies and shows up below 5 x 10"^ Hz. 
The subtraction of uniform baselines from the noise stream has 
a chopping effect that transforms a part of the low frequency, 
/ < 1/ W, power into high frequency, / > 1/ W- 
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Fig. 19. The effect of subtracting both the white noise and the 
1 // baselines on the power spectrum of the 1 // noise. The solid 
black line is the original 1 // spectrum. The solid colored lines 
show the power spectra of - FAn and the dashes lines the 
power spectra of - FRn. 



Fig. 20. Pixelization noise. In the top panel we show a piece 
of the original signal stream s (black) and the pixelized signal 
stream PBs (red). Their difference Zs = s - PBs is the pixeliza- 
tion noise {bottom panel). 



1 



2^base 

for which 



1 + 



1 



5 - a 



II a 



fk 



IJa 



and 

^2 



a 

cr^ \5 - a) a - 1 

For a = 1 this becomes 
1 



a-1 



2^ba 



= fk 



'-.mm. 

4 fklfc 



(58) 



(59) 



(60) 



(61) 



For fk = 50 mHz and a = l.l this gives ^base = 8.93 s, for 
which Vo^^To^ = 0.05936cr = 160.3yuK. This result does 
not take into account baseline errors and signal baselines (see 
Sect.EU). 

However, because of the different correlation properties of 
the different residuals, they have a different impact in map- 
making, and it is not enough to consider the time-domain vari- 
ance. So we need to return to this issue later, after we have stud- 
ied the residuals in the map domain and their angular power 
spectra. 



5. 7. Optimal baseline length 

What is the optimal baseline length to use? A partial answer can 
be found by minimizing the variance of the residual correlated 
noise, i.e., the residual noise minus the original white noise, g 

nc - FAn = (n^ - FRn^) - FRw - F(A - R)w - F(A - R)nc(57) ^ 

a; 

The variances of the four components are cr^, cr^^, cr^^, and cr^^. g 
The first three components are uncorrelated, so that their total g 
variance is cr^^ -h cr^^ -h cr^. Since cr^^ and cr^^ are much smaller 
than (T^^ -h cr^, we make the approximation where we drop the 
last two components. 

Using Eqs. (|49l) and ([53]) we find that cr^^ -h crl is minimized 

for 




0.00 



0.02 



0.04 0.06 

lag [s] 



0.08 



0.10 



Fig. 21. Autocorrelation function of the pixelization noise. Note 
the alternating correlation and anticorrelation for small lags, and 
the correlations at 1 minute lag, when the scanning returns to the 
same location on the sky. Due to spin rate variations the correla- 
tions are spread around the 1 minute value, and more around 2 
minute. 



5.8. Pixelization noise and signal baselines 

Pixelization noise, Zs = s - PBs, arises from signal gradients 
through a combination of pixelization, scanning strategy, and 
sampling frequency. These lead to correlations between close-by 
samples, and also correlations between samples from the same 
locations on the sky which are not close to each other in time 
domain. 

In Fig. [20] we show a short piece of the pixelization noise. 
Its autocorrelation function is shown in Fig. [21] We see that 
neighboring samples (lag = 1 = (1/76.8) s) are anticorrelated, 
whereas we have a positive correlation for lag = 2. The power 
spectrum is close to that of white noise, except that there are 
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some features near the Nyquist frequency due to these correla- 
tions. 

Comparing the pixel size Op = 6.87' to the sample separation 
6s = 4.68' and remembering that the scanning direction is mostly 
close to the direction of the pixel diagonal (both are often close 
to the ecliptic meridians), so that the pixel geometry tends to 
repeat at ~ = 9.1V intervals, close to 26>, = 9.38', we 

see that a pair of samples with lag 2 tends to land in about the 
same relative location within their respective pixels. With this 
small separation there is a positive correlation in the CMB signal 
gradient (smoothed with the beam). These combine to give a 
positive correlation between lag 2 samples. On the other hand, 
neighboring samples often land at opposite sides of the same 
pixel, or of neighboring pixels, leading to a negative correlation 
in their pixelization noise. These correlations would be different 
for different ratios of sample separation to pixel size. 

Since the pixelization noise arises from the signal, we ana- 
lyze it in the combinations s+ = (s^ -h S|,)/2 ("temperature") and 
s_ = {Sa - Sb)/2 ("polarization"), so that 



and 



s- = Q cos 2i^a + U sin 2i//a . 



(62) 



Each of the s+ and S- time streams contains nt/2 samples. For the 
ideal detectors considered here, s+ arises from gradients in the 
intensity of the signal, and s_ from gradients in the polarization 
of the signal. We assume the signal is statistically isotropic, and 
that Cf ^ = 0. 

We find that the expectation value for the variance of the 
pixelization noise is 



where 



rit 



p tep 



A : 



-E— E^'-- 



(63) 



(64) 



tt'ep 



and = - u^. Here is the unit vector giving the pointing 
(6t, (pt) of sample and is the unit vector pointing to the center 
of the pixel p hit by sample t. The approximation in Eq. ([63l) is 
the small-pixel approximation |r^| ^ 1. For (Tp(pol) we have 
also assumed an optimal distribution of polarization directions 
if/t within each pixel. 

This result, Eq. ([63]) , can be compared to the variance of the 
signal itself 

cr2(temp) ^ (/(n)/(n)) = ^ ?lllc[^ 



cr^(pol) ^ <e(n)G(n)) = <^(n)t/(n)) 



(65) 



Assuming the pixels are perfect squares, in the limit of a 
large number of hits ^hit,/? uniformly distributed over the pixel 
we get 




I 40.0 fjK 



4o:o 



Fig. 22. Pixelization noise rms in map pixels for the two 10° x 
10° regions. 



Since HEALPix pixels are not square, but somewhat elongated, 
we expect the actual F - A to be somewhat larger. In principle 
it can be calculated from the pointing data for a chosen pixeliza- 
tion. 

For our input Q, setting F - A = 0^/6, Eqs. ([63l) and ([65l) 
give cr,(temp) = 100.792yuK, crp(temp) = 13.332yuK, o-,(pol) = 
2.530yuK, and cr^(pol) = 0.781 yuK. 

The actual pixelization noise level in the simulation was 
(Tp(temp) = 13.516yuK and (Tp(pol) = 0.776yuK. 

By definition, a binned map of pixelization noise vanishes, 
BZs = Bs - BPBs = 0. Instead we can make a map of the rms of 
the pixelization noise at each pixel, by squaring each element of 
Zs, binning the resulting time stream into a map, and taking the 
square root of each / pixel. See Fig.[22l Comparing to Fig. [3] we 
see that pixels of large pixelization noise tend to "outline" hot 
and cold spots of the signal. 

Signal baselines arise from the pixelization noise in the same 
manner as white noise baselines arise from white noise, and 
we can divide them into the reference baselines of pixeliza- 
tion noise, and signal baseline errors. As = RZs -h (A - RZ)s. 
Approximating the pixelization noise as white, we get an esti- 
mate for the standard deviation (stdev) of RZs 

(Tn(temp) ..^ 
cr^,(temp) ^ ^ 1.521 yuKw(s)-'/' 



V^base 

cr^r(pol) ^ ^^^E^^ 0.0891 yuKw(s)-'/'. 

V^base 



(67) 



In Table [7] we show the stdev of the signal baselines As and 
their two contributions RZs and (A-RZ)s, and compare them to 
the estimate ([67l) for RZs. We see the estimate is quite good for 
RZs. The estimate is slightly off mainly because the actual pix- 
elization noise was slightly larger for temperature and slightly 
smaller for polarization than the analytical estimate (|63]) . 

In Fig. [23] we show the autocorrelation functions of the refer- 
ence baselines of pixelization noise and the signal baseline error. 
Unlike for white noise, now also the reference baselines are cor- 
related. This correlation extends over several repointing periods. 
These correlations enhance the relevance of the signal baselines 
in the map domain. 



6. Map domain 



and 



A«F 



Destriping is a linear process. The output map can therefore be 
(66) viewed as a sum of component maps, each component map be- 
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pixelization noise reference baselines 




time [h] 

pixelization noise baseline misestimation 




Fig. 23. Correlations in the signal baselines. The top panel 
shows the autocorrelation function of the reference baselines 
RZs of pixelization noise, and the bottom panel that of the sig- 
nal baseline error (A - RZ)s. The colors and line styles are as in 
Fig.im 



^base 


Eq. (|67l) 


RZs 


(A - RZ)s 


As 


As (G) 


2.5 s 


0.962 


0.972 


0.672 


1.297 




5s 


0.680 


0.686 


0.300 


0.831 




15 s 


0.393 


0.396 


0.121 


0.461 




1 min 


0.196 


0.200 


0.075 


0.248 


0.243 


1 h 




0.136 


0.075 


0.199 


0.196 


2.5 s 


0.0564 


0.0560 


0.0277 


0.0701 




5 s 


0.0399 


0.0396 


0.0145 


0.0476 




15 s 


0.0230 


0.0229 


0.0068 


0.0269 




1 min 


0.0115 


0.0115 


0.0043 


0.0144 


0.0147 


1 h 




0.0083 


0.0043 


0.0120 


0.0124 



Table 7. Statistics of the different contributions to the signal 
baselines. We give the standard deviation (square root of the 
variance) in jjK of the reference baselines of pixelization noise, 
the signal baseline error, and the solved signal baselines. The top 
half of the table is for the temperature stream {a + b)l2 and the 
bottom half for the polarization stream (a - b)l2. The column 
marked "(G)" was obtained using galactic coordinates for the 
HEALPix maps. 



0.01 




full, I-Stokes 
cmb 

white noise 
1/f noise 
full-components 
full, Q-Stokes 
full, U-Stokes 



100 

baseline length [s] 

Fig. 24. Linearity of destriping. The three thin black lines at the 
bottom shows the rms difference between the "full" map, ob- 
tained from the sum of the component TODs, and the sum of 
the component maps. This is compared to the signal and noise 
levels in the maps; different colors showing the different compo- 
nents and the thick black line the full map; different line styles 
correspond to the three Stokes parameters. 



ing a result of destriping one individual TOD component: signal, 
1 // noise and white noise: 



mo 



= B(s - FAs) + B(nc - FAn^) + B(w - FAw) . 



(68) 



As the three TOD components are statistically independent, 
so are the corresponding component maps. The expectation 
value of the map rms is thus obtained as the root sum square (rss) 
of the expectation values of the rms of the component maps. 

In this paper we have considered a single 1-year realization 
of each TOD component only. There are likely to be random 
correlations between the component maps. However, the output 
map and residual map rms we get when we produce maps from 
the full TOD, agree to better than 0.5% with the rss of the rms 
of the corresponding maps produced from component TODs. 

6.1. Linearity 

In practice the linearity of destriping is affected by the numer- 
ical accuracy of the baseline calculation, which depends on the 
convergence criterion for the conjugate gradient iteration. We 
checked the linearity by comparing the sum of the component 
maps to a map obtained directly from the sum of the compo- 
nent TODs. The difference is well below the nK level, except for 
baselines shorter than 10 s. See Fig.[24l 

6.2. Splitting the residual map into components 

The three component maps can each be further divided into a 
binned map and a baseline map. 

The correlation properties of these components are different 
in the map domain from the TOD domain. 

While the white noise baselines FAw are correlated with the 
white noise w in the TOD, the correlation vanishes in the map: 
Writing 

Bw = M-^P^C;^w (69) 
BFAw = M-^P^C;^FD-^F^C;^ (l - PM-^P^C;^) w (70) 
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Map residuals 
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Fig. 25. The residual temperature map rms as a function of base- 
line length. The top panel is for niout -Bs, which is dominated by 
the binned white noise map, whose rms is 55.19 yuK independent 
of ^base- The bottom panel shows the other three components, 
and their root sum square, which is the rms of the residual map 
niout - Bs - Bw. 



10 




10000 



baseline length [s] 

Fig. 26. Same as the bottom panel of Fig. [25l but now we have 
further separated the residual 1// noise into unmodeled 1// 
noise {dashed) and 1// baseline error (dash-dotted)', the white 
noise baselines into reference baselines (dashed) and baseline 
error (dash-dotted); and signal baselines (shown here in green) 
into reference baselines of pixelization noise (dashed) and base- 
line error (dash- dotted). The dotted lines show the rss of the 
dashed and dash-dotted lines. For white noise it falls on the solid 
line, showing that the two white noise residual map components 
are uncorrelated. For 1 // noise and pixelization noise the two 
components are positively correlated. The black dotted line is 
the analytical approximation Eq. ([86l) . 

we get that 

<(Bw)(BFAw)^) = M-^P^C;\ww^)(l-C;^PM-^P^) 

C;^FD-iF^C;^PM-i . (71) 
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baseline length [s] baseline length [s] 

Fig. 27. The rms of the residual /, Q, and U maps, and their 
different components. The overall residual power in a destriped 
polarization map is higher than in a temperature map. The CMB 
signal residual in turn is much weaker. Note how Q is larger 
than U in the signal baselines. Smoothing the residual maps with 
a Gaussian beam (dashed lines) removes excess power at sub- 
beam scales. 




-i2Q ^^^^^ 120 p.K 



Fig. 28. The residual / map including binned white noise, niout - 
Bs, for ^base = 1 min. It looks the same for other baseline lengths. 
Also the binned white noise map Bw looks the same, as it is the 
dominant component in the residual map. The Q and U maps 
look the same, but have a larger amplitude. 



By substituting <ww^) = one readily sees that the correlation 
vanishes, 

((Bw)(BFAw)^) = . (72) 

Because the binned white noise map and the white noise 
baseline map are independent, the (expectation value of) the rms 
of the residual white noise map Bw - BFAw is obtained as the 
rss of the rms of the two maps. 

This means, that while subtracting the white noise baselines 
removes power from the TOD, it adds power to the map. If the 
noise is pure white, naive binning produces a better map than 
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destriping. On the other hand, baselines of the 1 // noise are cor- 
related with the 1 // noise itself both in the TOD and in the map. 

We can use the rms value taken over all the residual map 
pixels as a figure-of-merit for the map-making method. We cal- 
culate it separately for the three Stokes parameters, / (tempera- 
ture), 2, and U. Note that we subtract the I monopole from the 
residual map before calculating the rms, since it is irrelevant for 
CMB anisotropy and polarization studies, and destriping leaves 
a spurious / monopole in the map. See Figs.[25l[26l andlTTl 

The difference between the output map and the binned signal 
map 



- Bs = Bw - BFAw + B (lie - FAn^) - BFAs 



(73) 



is the residual map including binned white noise. We see that it 
can be divided into four uncorrelated contributions: the binned 
white noise map, the white noise baseline map, the residual 1 // 
noise map, and the signal baseline map. 

The niout - Bs / map is shown in Fig.[28l This map is dom- 
inated by the binned white noise map, which is independent of 
the baseline length. Therefore the visual appearance of the resid- 
ual map is the same for all baseline lengths, as it looks the same 
as the white noise map. Also the residual Q and U maps look the 
same, just with ~ ^^2 larger amplitude. 

The binned white noise map is independent of the baseline 
length, and it is an unavoidable component of the residual map, 
uncorrelated with the other components. Therefore we focus on 
the rest of the residual map, without the binned white noise com- 
ponent, i.e.. 



mo 



Bs - Bw = -BFAw + B (lie - FAn^) - BFAs . 



(74) 



For the rest of this paper, the residual map without further qual- 
ification, refers to this map. See Figs.[29land[3Ql 
We divide this further into 



mout - Bs - Bw = -BFRw - BF(A - R)w 

+B(iie - Fa,) - BF(Aiie - a,) 
-BFRZs - BF(A - RZ)s . 



(75) 



in analogy with Eq. (|48l) . Of these six components, the unmod- 
eled 1 // map B(iic - Fa,) is correlated with the 1 // baseline er- 
ror map BF(Aiic - a,), and the pixelization noise reference base- 
line map BFRZs is correlated with the signal baseline error map 
BF(A - RZ)s. Otherwise the components are uncorrelated with 
each other. The division is useful, since the diff'erent components 
have a diff'erent structure in the map domain. The unmodeled 1 // 
map consists of diff'erent noise frequencies with half- wavelength 
mostly of the same order or shorter than the baseline (Fig. [15]). 
The other 5 components consist of baselines laid on the map. 
For the white noise reference baseline map, these baselines are 
uncorrelated with each other, the other 4 components have dif- 
ferent levels of correlation. In Fig.[3T]we show each of the six 
components for a 10° x 10° region of the / map. 

In Fig.[26l we show separately the rms of these components. 

6.3. Analytical estimates 
6.3.1 . White noise baselines 

Since the white noise reference baselines are uncorrelated with 
each other, the white noise reference baseline map rms can be 
estimated. If all hits to a pixel came from a diff'erent baseline, 
we could treat them as white noise, with the white noise refer- 
ence baseline variance cr^ = cr^/^base- Thus the variances of !„, 
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Fig. 29. The residual / map (with binned white noise sub- 
tracted), niout - Bs - Bw, for diff'erent baseline lengths (from 
top down:) 1 hour, 1 minute, 15 seconds, and 2.5 seconds. The 
Q and U maps look the same with ~ ^/2 larger amplitude. 



Qp, and Up would be just (Mp)-\l, l)/^base ~ Cr^/^base/^hit,p, 
(Mp)-H2,2)/^base ~ 2o-V^baseMhit,/;, and (Mp)-H3, 3)Mbase ~ 
2(T^/^base/%it,/? . 

For baseline lengths less than or equal to the spin period, 
^base ^ Imin, this holds if the sample separation is larger than 
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Fig. 30. Same as Fig. [29] but for the 10° x 10° region near the 
ecliptic. 



the pixels, 6s » Op. However, in our case 6s < Op, and two or 
three successive samples may hit the same pixel. These succes- 
sive hits are then almost always from the same baseline, and, for 
the baseline components, fully correlated, i.e., equal. 

Denote by nfn the fraction of hits to a pixel that belong to a 
sequence of exactly n successive hits to the same pixel, i.e, there 
are /«^hit,/? such sequences, and Yjn ^fn = 1- The variance of Ip 
is then 



{Il} = \j](n'fn) 



(76) 



where 2:„(nV«) > 1- 

With Eq. $4% we have for the I map rms 







\ nut 1 


V^base 



(77) 



To get an estimate for (Y^ni^^fn)) we assume a square pixel 
and a uniform hit probability distribution within the pixel area. 

For scanning in the direction of the pixel diagonal, for 6s > 
^^26p there are no successive hits to the same pixel and /i = 
Tuni^^fn) = 1- For ^p/ ^ <0s < ^Op there can be a maximum 
of two hits, with (/2) = (l - r/ ^/2f and 



(78) 



and with V2^p/3 < 6s < Op/ ^/2 3. maximum of three hits, with 
</i) = r2,(/2) = 3V2r-l-7r2/2and 



^(^2^„) = 5-6V2r + 5r^ 

n 

where r = OJOp. 



(79) 
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Fig. 31. The six components of the residual / map shown for 
the 10° X 10° region near the ecliptic (^base = 1 min). Top left: 
Unmodeled l/f noise B(nc - Far). Top right: l/f baseline er- 
ror BF(Anc - sir). Middle left: White noise reference baselines 
-BFRw. Middle right: White noise baseline error -BF(A-R)w. 
Bottom left: Pixelization noise reference baselines -BFRZs. 
Bottom right: Signal baseline error -BF(A - RZ)s. 



For scanning in the direction of the pixel side, two hits, but 
no more, are possible for Op/ 2 < Os < Op with /2 = 1 - r. 

For our case, V2/3 < 1 /2 < r = 0.682 < 1 / V2, which gives 
^(^V.)) = 1.54 = (1.24)2 



(80) 



for scanning in the diagonal direction, and {Yjni^^fn)) = 1.64 = 
(1.28)2 for the pixel side direction. 

However, HEALPix pixels are not square, but can be signif- 
icantly elongated. In principle, the mean value ( ^n^^" ) could be 
calculated from the pointing data and the chosen pixelization. 
Here we just take it by comparing the actual /rms from the maps 
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^base 


-^rms 






^rms 


^base 


-^rms 


Qrms 


^rms 


2.5 s 


4.947 


1.542 


6.974 


7.072 


2.5 s 


6.729 


9.303 


9.580 


15 s 


2.006 


1.521 


2.846 


2.914 


15 s 


2.319 


3.291 


3.356 


1 min 


1.003 


1.522 


1.435 


1.466 


1 min 


1.181 


1.678 


1.715 


1 h 


0.895 




1.289 


1.320 


Ih 


1.088 


1.552 


1.590 



Table 8. Statistics of the white noise reference baseline maps. 
The map rms values given are in yuK. The third column is an 

estimate of ( ^^^^/" )/(^i^t) obtained by comparing I^ms to Eq. (TTTI) . 



^base 


-^rms 


c 


Qrms 


^rms 


2.5 s 


4.543 


0.168 


6.223 


6.433 


15 s 


1.165 


0.092 


1.664 


1.666 


1 min 


0.634 


0.102 


0.881 


0.904 


1 h 


0.630 




0.875 


0.898 



Table 9. Statistics of the white noise baseline error maps. The 
map rms values given are in jjK. The third column is an es- 
timate of the correlation between baselines contributing to the 
same pixel obtained by comparing /rms to Eq. ([83]) . 



Table 10. Statistics of the white noise baseline maps. The map 
rms values given are in ^K. 

^base Eqs. (1561841) Eq. ([84]) I,^, Q^, U^, 

7.5 s 2.994 3.081 3.062 4.314 4.358 

15 s 3.727 3.901 3.852 5.450 5.526 

Imin 6.053 6.313 5.778 8.128 8.225 

Ih 9.357 13.02 13.178 

Table 11. Statistics of the unmodeled 1 // noise maps. The map 
rms values given are in jjK. The second and third columns are 

estimates for the / rms based on Eq. (1841) , with ^ ^nj'' ) taken 

from Table m The third column uses the actual cTu, whereas the 
second column uses the analytical estimate (l56b . 



(see TablelSl) to Eq. dTTl) . This gives (from ^base = 15 s and 1 min) 
(Z^VJ^hit) ~ 0.000636 ^ 1.52<^-J) ^ (0.025)^. We denote 



"hit 



Wt> - 1.52. 



(81) 



Eq. dSOb is quite close to this. 

For the other Stokes parameters, we expect 



^rms ~ <^rms 



V2/n 



(82) 



where the approximation corresponds to assuming an ideal dis- 
tribution of polarization directions. We expect this approxima- 
tion to be good for the full year data, since most pixels have 
rcond values close to 0.5. 

For baselines longer than the spin period, contributions to a 
pixel from successive scan circles tend to come from the same 
baseline, so the reduction in the white noise baseline variance 
is canceled by the reduction in the number of contributing base- 
lines. Thus the map rms from white noise baselines is almost flat 
between ^base = 1 min and ^base = 1 h. 

Although in the time domain the white noise baseline error is 
much smaller than white noise reference baselines, their corre- 
lations make them important in the map domain. Assuming the 
correlation between (a -h b)l2 baselines contributing to the same 
pixel were c, the expected variance of the white noise baseline 
error map would be 



2cr 



1 -h 



(83) 



where o-^e+ is the stdev of the {a + b)/2 white noise baseline 
errors, given in TableO Assuming Eq. ([83l) to hold for our maps, 
we have solved for c for difl'erent ^base in Table [9l These numbers 
can be compared to Fig. [TT] and Tabled 

Since the white noise reference baselines and the white 
noise baseline errors are uncorrected, the full white noise base- 
line map rms is close to the rss of these two components. See 
Tablefini 



6.3.2. Unmodeled 1// noise 

Most of the power in unmodeled 1 // noise is in frequencies near 
l/(2^base)- Therefore, for 6p «: ^base successive hits to the same 
pixel should be almost fully correlated. When ^base is much below 
the spin period, hits from difl'erent spin periods should be almost 
uncorrelated. Thus we can estimate the unmodeled l/f noise 
map rms in the same manner as the white noise baseline map, as 



\ nut 



(Tu ~ 0.025cr^ . 



(84) 



When ^base is comparable to the spin period, there will be cor- 
relations (positive or negative) between hits from nearby spin 
periods. We compare Eq. ([84l) to the actual binned maps of un- 
modeled 1 // noise in Table [TT] We see that for ^base = 1 min, 
Eq. ([84b is an overestimate, indicating that there are negative 
correlations between hits from nearby spin periods. 

6.3.3. Ideal scanning 

We define ideal scanning so that the pointings from the diff'er- 
ent scan circles of the same repointing period fall on top of each 
other, i.e., there is no nutation and the sampling is synchronized 
with the spin period. In this case, the part of the unmodeled 1 // 
noise for long baselines (^base a multiple of the spin period) that is 
modeled by 1 min baselines gets totally averaged out, so that the 
contribution from unmodeled 1 // noise to residual maps would 
stay constant from ^base = 1 min to ^base = 1 h. In our nonideal 
case, some of this noise leaks out, so that the unmodeled con- 
tribution rises slowly in this range also. See dashed red line in 
Fig. [261 

Likewise, for an ideal scanning, the white noise reference 
baselines make an equal contribution to the map for any baseline 
length that is an integer multiple of the spin period, and fits into 
the repointing period an integer number of times. 

6.3.4. Total noise 

From Fig. [26] we see that the two dominant contributions to the 
residual maps are the white noise reference baselines and the 
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unmodeled 1 // noise. For both of them we have analytical esti- 
mates, and both of them map from time domain to map domain 
in roughly the same way. Thus we get an analytical estimate for 
the residual / map rms by multiplying the 



estimate 

from Eqs. (|49]) and (|53 with ^/(^^) ~ 0.025 for ^base < 1 

min. When > 1 and ^base » Sample, so that {fxifcf ^ 1, we 
have 



1 



\ /sample 



/sample ^base 

4 



{5-a){a-\) 



(85) 



This gives for our case {cr = 2700yuK, ((^ n^fn)} = 1.52, n'J = 
0.000418, /sample = 76.8 Hz, a = 1.7), 



I' 



(86) 



where ^base is given in seconds. (For ^base ^ 1 ™n, our analytical 
treatment can just estimate that the map rms should stay constant 

from ^base = 1 min to ^base = 1 h.) 

However, this estimate is not as good as in the time domain, 
since the importance of the neglected components, the baseline 
errors, has grown dramatically when going from the time domain 
to the map domain. See Fig. [26l Since these components rise 
towards shorter baselines, the residual map rms is minimized at 
a somewhat larger ^base than Eq. (|58]) gives. 

6.4. Pixelization noise and signal baselines 

For pixelization noise, already the reference baselines are 
strongly correlated (see Fig. [55]) and therefore their map rms 
cannot be estimated like for white noise reference baselines and 
unmodeled 1 // noise. Due to these correlations their impact in 
the map level is significantly larger than their small variance in 
the time level (see Table [7]) would indicate. Instead, for both the 
pixelization noise reference baselines and the signal baseline er- 
ror, the situation is similar to white noise and 1 // noise baseline 
error. 

For the residual l/f and white noise baselines, the Q and 
U maps look the same as the / maps, just with a factor ~ 
V2 larger amplitude, since they originate from the same time- 
domain noise, which is independent for each detector. 

For the signal baselines (see Fig. [32l) the situation is, how- 
ever, diff'erent, since they originate from the signal, where Q and 
U are much smaller than /. We also note that Q is much larger 
than U, although they are of same magnitude in the signal. This 
is related to the coordinate dependence of the definition of the 
Stokes parameters Q and U together with a combination of fac- 
tors in our study. First, the signal contains only E mode polar- 
ization, which means that Q has structures along the coordinate 
lines, whereas U has structures oriented 45° from them. Second, 
we are using ecliptic coordinates, and we have employed a scan- 
ning strategy, where the scanning goes almost parallel to the 
lines of longitude for a large part of the sky. The signal baselines 
originate from the signal gradients within pixels. For a signal 
structure oriented along the scanning direction, the signal gradi- 
ent structure remains similar for a sequence of pixels along the 
scanning. Thus the measurement diff'erences between diff'erent 
scans through these pixels are similar for a sequence of pixels, 
favoring their misinterpretation as noise baselines. 

To verify the eff'ect of the coordinate system, we redid the 
^base = 1 h and 1 min cases using galactic coordinates. See 
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Fig. 32. The signal baseline / (top), Q (middle), and U (bottom) 
maps -BFAs (^base = 1 min) for the two 10° x 10° regions. 



coord. / rms Q rms U rms P rms 



1 min 


E 


123.3 


9.4 


4.5 


10.4 


1 min 


G 


115.9 


6.8 


7.9 


10.4 


Ih 


E 


124.0 


9.4 


4.6 


10.5 


Ih 


G 


116.5 


6.9 


7.9 


10.5 



Table 12. Eff'ect of the coordinate system (E = ecliptic, G = 
galactic) on the signal baseline maps. The map rms is given 
in nK (not in /iK, like the other tables!), and P stands for 



Table [121 We see that the asymmetry between Q and U largely 
disappears, but the total polarization signal residual is not much 
aff'ected. For the temperature residual we see a small improve- 
ment. This is partly explained by the reduction of the signal base- 
line variance, seen in Table [71 

6.5. Residuals at different angular scales 

The residual map rms alone is a poor measure of the quality of 
the output map. Since the nature of the residual (see Figs. [29l and 
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[30l) is different for different baseline lengths, we need to look at 
the structure of the different map residuals in more detail. 

For long baselines, the residual mostly comes from the part 
of the 1// noise that cannot be modeled with baselines, and 
appears mostly at very small angular scales on the map, near 
the pixel scale; whereas for shorter baselines it comes from un- 
wanted baselines, which appear as larger scale structures. This 
can be seen from Fig.[27l where we have smoothed the residual 
map with the detector beam, before taking the rms. This smooth- 
ing almost erases the difference between the 1 // residuals for 
baseline lengths from ^base = 1 min to 1 hour. This is because the 
1 min scanning circles fall almost on top of each other during the 
1 hour repointing period, and the width (due to nutation) of the 
ring on the sky traced by the beam center during a repointing pe- 
riod is less than the beam width. Beam- smoothing has much less 
effect on the white noise baseline map and the signal baseline 
map. Baseline lengths ^base = 10s to 15s still give the smallest 
total residuals, but the difference from longer baselines is much 
reduced by beam- smoothing. Since residuals at larger scales are 
for most purposes more harmful than sub-beam residuals on the 
map, we conclude that it is better to choose a somewhat longer 
baseline than what would minimize the residual map rms. 

6.6. Angular power spectra of map residuals 

Since we are considering full-sky maps, their angular power 
spectra Cf^ can be calculated directly from them (we used 
anaf ast of the HEALPix package). 

We plot the angular power spectra of the residual maps in 
Figs.[33]and[34]for different baseline lengths. It is clear that base- 
lines shorter than ^base = 10 s, lead to more large scale structure 
in the residuals. Long baselines lead to a high-^ tail in the resid- 
ual that appears much like white noise (flat Q). 

In Figs. [35] and [36l we show angular power spectra Cj^ and 
Cf^ of different map components for the cases ^base = 1 h, Imin, 
15s, and 5 s. 

For baselines that are multiples of the spin period (1 min 
and 1 h) we see the characteristic even-odd variation in the Q 
of the baseline components. If the scanning circles had the full 
90° radius, they would contribute only to the even multipoles. 
In our case the circle radius is 87.77° = 90° - 2.23°, and there- 
fore we see a beat pattern, where the maximum even-odd multi- 
pole difference occurs at multipoles f that are near multiples of 
90°/2.23° = 40.4. For the low £ of the unmodeled 1 // contribu- 
tion we see the opposite pattern, since the unmodeled 1 // con- 
tains mostly frequencies which vary just over those timescales 
over which the baseline contributions are constant. 

The angular power spectrum of full-circle uncorrected base- 
lines goes as 



TT 



.00 



Cf oc 



1 



2^-h 1 



(87) 



(Eftstathiou 12005b . Therefore we have plotted (2£ -h l)Q/7r in 
Figs. [35] and [36] We see that Eq. dSTl) indeed holds well for 
white noise reference baselines; even for ^base ^ l//spin, al- 
though these have less power at the lowest £. Although the un- 
modeled 1 // contribution does not consist of constant baselines, 
the correlations of the parts between different baseline segments 
are weak (nonexistent for ^base ^ l//spin). Therefore Eq. (ISTl) 
holds fairly well for the unmodeled 1 // also; except at low £ for 
short baselines, where there is a lack of power since the unmod- 
eled 1 // varies more rapidly along the scan path; and for high £ 
for ^base » 1 //spin, which havc excess power at high related to 
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Fig. 33. Angular power spectra of the residual temperature maps 
for different baseline lengths: 1 h (black), 4 min (blue), 1 min 
(light blue), 30 s (light green), 15 s (green), 7.5 s (pink), and 
5 s (red). The black horizontal line is the white noise level 
and the black dashed line is the theoretical CMB input spec- 
trum smoothed with the beam and pixel window functions. Top: 
Residual map including white noise, niout-Bs. Bottom: Residual 
map niout - Bs - Bw (binned white noise subtracted). 



imperfect superposition of the different scan circles of the same 
ring, mainly due to nutation. 

The other contributions have different angular scale depen- 
dencies, related to the correlations between baselines. We see 
that the baseline error components have steeper spectra than the 
reference baseline and unmodeled 1 // contributions. This makes 
them important at large scales (low multipoles), where they are 
comparable or even stronger than the white noise reference base- 
line and unmodeled 1 // components, which dominate at high £ 
and contribute most to the residual map rms. If one considers 
just the signal baselines, the baseline error completely dominates 
over the reference baselines for short ^base and low £. 

For long baselines (~ 1 min or longer), the unmodeled 1 // 
noise dominates the residuals for ^ > 10, but the l/f baseline 
error contribution can be comparable for ^ < 10. For shorter 
baselines, the white noise baselines become more important. 

The C^^ spectra of different residual components look qual- 
itatively like Cf^, except for the signal baseline components. 
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Fig. 34. Same as Fig. [33l but for the polarization E and B mode 
spectra. 



which have less power, reflecting the lack of B-mode signal in 
the input. Therefore we have not plotted the C^^ spectra, except 
for these signal baseline components, which we have included 



in Fig. [36] along with the Cf^ spectra. 

Because of large ^-to-^ variations in the residual Q, these 
plots are difficult to read. Therefore we also plot (Fig. [37]) square 
roots of the cumulative angular power spectra. 



2^+1 
An ' 



IV + 1 



{Cf ^Cf), (88) 



which give the total contribution to the residual / and F map rms 
from multipoles up to and 



\ 



V2r + l|^,^l 



471 



(89) 



The beam fwhm 6 = 12.68' corresponds to multipole £ = 
180°/^ ~ 850, so we are more concerned about the behavior 
up to this f than above it. 

Fig.[37]provides probably the most concise meaningful com- 
parison of the quality of maps vs baseline length. The ^base = 15 
s case appears the best in terms of cumulative residual power in 
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Fig. 35. Angular temperature power spectra of diff'er- 
ent map components for four different baseline lengths: 1 h 
{top panel), 1 min, 15 s, and 5 s {bottom panel). Solid black: 
Residual map (with binned white noise subtracted). This can 
be split into the following six components: Red: Unmodeled 
1 // noise. Pink: Additional effect of 1 // baseline error. Blue: 
White noise reference baselines Light blue: White noise base- 
line error Green: Reference baselines of pixelization noise. Light 
green: Additional effect of signal baseline error. The smooth 
black curve is the white noise level. 
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Fig. 36. Same as Fig.[35]but for the E mode polarization spec- 
trum Cf^. We also show the two signal baseline contributions of 
the B mode polarization spectrum C^^. They are the lower blue 
curves: Blue: Reference baselines of pixelization noise. Light 
blue: Additional effect of signal baseline error. 
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Fig. 37. Square roots of the cumulative angular power spectra 
(see Eqs.[88]and[89]) for the residual temperature (top panel) and 
polarization (middle) maps and TE correlation (bottom), up to a 
given multipole plotted as a function of £ for different baseline 
lengths: 1 h (black), 4 min (blue), 1 min (light blue), 30 s (light 
green), 15 s (green), 10 s (green dashed), 7.5 s (pink), and 5 s 
(red). 



the map for the relevant multipoles. Although ^base = 10 s pro- 
duces a smaller residual map rms, it is only because of its small 
sub-beam- scale residuals. Interestingly, the ^base = 1 h base- 
line seems to be the best for minimizing residual temperature- 
polarization correlations at intermediate scales. It is also better 
overall than the ^base = 1 min and 4 min cases for ^ < 300. 
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If one is only interested in large-scale features (low £) there 
are no big differences between any of the baseline lengths from 
15 s to 1 h, but 10 s or less should be avoided. For the lowest mul- 
tipoles there is some randomness in the results, since we studied 
only one noise realization, so one should not try to draw conclu- 
sions from the small differences seen there for ^base = 15 s to 1 h. 
(The noise residuals, for the cases ^base = 1 min and ^base = 1.25 s 
with noise prior, have been studied via Monte Carlo in Keskitalo 
etal. 120091 ) 

The signal baseline contributions appear a minor effect at 
all scales. In this study the signal contained only the CMB. 
In reality, the gradients in the signal are often dominated by 
foregrounds, and therefore the signal baseline effect is larger. 
Foreground signals are considered in Keihanen et al. (120091) . 
Foregrounds were also included in the map-making studies of 
Ashdown et al. (2007b, 2009 ), and especially in the former there 
was a detailed study on the signal baseline contribution and how 
it could be minimized. 

Map residuals influence the precision at which we are able 
to determine the angular power spectrum of the CMB map. 
We can subtract the expectation value of the Q of the resid- 
ual from the map spectrum, but individual realizations deviate 
from this expectation value, leading to an error in the CMB Q 
estimate. A multipole ^ of a map has at the best 2^ -h 1 degrees 
of freedom. Statistically isotropic signals (e.g. CMB or a white 
noise map of uniform pixel variance) have these degrees of free- 
dom. Deviations from the statistical isotropy may lead to corre- 
lations in the m-modes of a multipole which in turn decreases 
the degrees of freedom and therefore increase the error in Q. 
Fortunately the strongest correlations of the map residuals occur 
nearly along the ecliptic meridians. Correlations along meridi- 
ans do not lead to correlations in the m-modes. Therefore we 
can expect that the m-mode couplings due to map residuals are 
weak and the excess spectrum error is small. For now, we did not 
investigate these errors any further, but decided to leave this for 
future studies. 

6. 7. Low multipoles of I, Q, and U maps 

For cosmological purposes, one calculates the Cj^, Cj^, Cf^, 
and C^^ angular power spectra of the output maps, which repre- 
sent the fundamental properties of the temperature and polariza- 
tion field, and are coordinate independent. However, for analyz- 
ing residual map structure, it may be more intuitive to consider 
the Q and U maps as two separate maps of a scalar quantity, and 
calculate their ordinary (spin-0) angular power spectra. 

The Q and U are given in terms of the a^^ and a^^ as 
(Zaldarriaga & Seliak fTWTb 

Q^iU = -Y^{al^ial)2Y^, (90) 
where the 



V2f ■+■ 1 
-^e'""^di^_2(0) (91) 

are the spin-2 spherical harmonics (Newman & Penrose 1 19661 
Goldberg et al. 1967 ) and the ^, (0) are the Wigner J-functions 
(we follow Varshalovich et al. 1988 ), which are real. 

Thus the monopoles Q and U of the Q and U maps are given 

by 

Q + iU = -J]{al + ial)2n, (92) 




Multipole, I 

Fig. 38. The effect of destriping on the low multipoles of 
/ (blue), Q (red), and U (green) maps. The solid lines give 
the multipoles of the binned signal (CMB) map, where the 
monopoles of / and U, and the dipole of / should ideally vanish 
in our case. The dashed lines are for the naive binned map in- 
cluding the noise, and the dash-dotted lines are for the destriped 
(^base = 1 min) map. Destriping is not able to remove the noise 
monopole of the / map (the monopole does change but does not 
become small), but the noise monopoles of Q and U maps are 
removed, about equally well as the other low multipoles. The 
dotted lines are for the binned white noise map. The horizontal 
black lines give the expected white noise levels. 



where 

^Y^ ^ ^ f 2YT(0,cp)dQ 



d{_2(e)^mede, (93) 

which are real. Thus the Q and U monopoles are given by 

G = -Xi"^o2?? and C7 = -J]<2}^, (94) 

Since in our case the input spectrum contained no B mode, the 
monopole vanishes in the input U map. 

Fig. [38l illustrates the effect of destriping on the /, Q, and U 
maps. 

The low multipoles of the residual maps are shown in 
Fig. [39l divided by the white noise level. We see that the resid- 
uals at lowest multipoles are larger than the white noise level. 
The baseline length does not make a large difference for these 
low multipoles, except that too short baselines (less than 15 s, 
not included in Fig. [39]) should be avoided. 



7. Effect of noise knee frequency 

The importance of the correlated 1 // noise depends on its am- 
plitude and spectrum. In this paper we do not consider the effect 
of possible spectral features in the noise, and we have parameter- 
ized the noise just by the 1 // slope -a, knee frequency and 
white noise level. Since we produced the simulated 1 // noise 
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Fig. 39. The / , Q, and U low multipoles of the residual maps for 
different baseline lengths: 1 h (black, 4 min (blue), 1 min (light 
blue), 30 s (light green), and 15 s (dark green). For this plot they 
are divided by the white noise level. The top panel shows /, the 
bottom panel Q (solid) and U (dashed). 



Residual noise in the output map 



signal baselines 
white noise baselines 
residual 1/f (fj^=100mHz) 
residual 1/f (fj^=50mHz) 
residual 1/f (f =25mHz) 




Binned 1/f, f =100mHz 



Binned 1/f, f, =50mHz 



Binned 1/f, f, =25mHz 



Binned white noise level 



100 

baseline length [s] 

Fig. 40. The knee frequency directly relates to the amount of 
residual 1 // noise in the output map. Reducing the amount of 
correlated noise (and hence the knee frequency) causes the white 
noise baseline contribution to exceed residual 1 // contribution 
at longer baselines. Therefore lower fk means longer optimal 
baseline length. 
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Multipole, t 



Fig. 41. Same as Fig.[39l but for a lower knee frequency fk = 25 
mHz (and a smaller selection of baseline lengths). 



separately from the white noise, we can change fk simply by 
multiplying the 1 // part by 



n^2 



fk,2 
fk,l 



a/2 



(95) 



This will change the residual 1/f contribution to the residual 
map by the same factor, while the white noise baseline and sig- 
nal baseline contributions are unaffected (unless one changes to 
a different ^base). For angular power spectra the 1// contribu- 
tions are rescaled by the square of this factor. See Fig.|40l where 
we have plotted the residual map rms as a function of ^base for 
fk = 100 mHz and 25 mHz, besides the case of fk = 50 mHz, 
which we have so far considered. A higher fk favors shorter base- 
lines, since the stronger 1 // noise at relatively high frequencies 
needs to be modeled better. A lower fk, on the other hand, favors 
longer baselines, since then the white noise and signal baselines 
are relatively more important. For ~ 10 mHz the minimum 
rms would move to ^base ~ 1 min. 

In Fig.|40lwe show also the rms of the binned 1 // noise map. 
As we lower fk it moves down. For very small fk it would fall 
below the white noise baseline rms. At this point simple binning 
would produce a better result than destriping. For our simulated 
1 // noise, and for long baselines, this would happen at the ex- 
tremely low fk = 21.6yuK. It should be noted however, that in 
our case the binned 1 // rms is heavily dominated by the lowest 
frequencies, and in other cases (smaller slope a, larger /mm) the 
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Fig. 42. Same as the bottom panel of Fig. [25l but the solid lines 
correspond to using a noise prior. The dashed lines are from 
Fig. [25l The noise prior has practically no effect for ^base ^ 1 
min. 




Multipole, I 

Fig. 43. Same as top panel of Fig.[37]but for the case of using the 
noise prior. The colors and linestyles are the same as in Fig.[37l 
except we have added the ^base = 2.5 s case (solid red), and the 
^base = 5 s casc appears now as the dashed red line. 



relevant 1// rms could be smaller for a given so that simple 
binning could become superior already at a higher 

In Fig. Si] we show the low multipoles of the residual /, 2, 
and U maps recalculated for fk = 25 mHz. 

8. Effect of noise prior 

Although short baselines can potentially model correlated noise 
better, they fail because of the large random amplitudes they pick 
from white noise. This can be remedied by applying prior infor- 
mation on the noise spectrum to prevent too large differences 
between amplitudes of nearby baselines. This is discussed in de- 
tail in Keihanen at al. ( 2009) , but we give a short preview of the 
results here. See Figs. |42l |43l and|44l The noise prior has little 
effect for ^base » Vfk, but for short baselines the effect is dra- 
matic. We note that now the results keep improving as the base- 
line is shortened, at least until ^base = 2.5 s, the shortest we tried. 
(Using very short baselines with the noise prior makes the code 




2 4 6 8 10 

Multipole, t 



Fig. 44. Same as Fig. [39l but with the noise prior, and using 
a short baseline ^base = 2.5 s. The Stokes parameters are now 
indicated by the colors: blue for /, red for Q, and green for U. 
The solid lines are for fk = 50 mHz, and the dashed lines for 
fk = 25 mHz. 



I map rms 

1 5 I ' ^V, ' ^ 




1 10 100 1000 10000 

baseline length [s] 

Fig. 45. Same as the bottom panel of Fig. [25l but now we have 
included results from 7-month (dashed) and 1 -month (dash- 
dotted) survey segments also. The 1 -month results show aver- 
ages from all 12 months of the survey. 



more resource intensive.) For the very lowest multipoles, the re- 
sults with the short baselines do not, however, become much bet- 
ter than the ones obtained with longer baselines (with or without 
noise filter). Compare Fig.|44lto Figs. [39l and |4T] 

9. Maps from shorter survey segments 

We now consider making maps from a shorter time segment of 
the data. In a full year of observations, for a PLANCK-like scan- 
ning strategy, all parts of the sky are looked at two different sea- 
sons. About 7 months is needed to observe the full sky. Maps 
from a shorter segment cover just a part of the sky. 

When maps are made from a shorter observation period, the 
number of crossing points is reduced and their pattern is differ- 
ent. Destriping 7 months of data differs qualitatively from the 
case of a full year, since for a large part of the sky the second- 
season observations are missing. We expect the loss of the cor- 
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Fig. 46. Linear hit maps near the ecHptic North Pole for the 
one-month surveys corresponding to the first (left) and second 
(right) month of the full year survey. Note how diff'erent is the 
structure of crossing points between these two months. The same 
happens near the South Pole. In the first month all circles cross 
each other at almost the same pixel, whereas in the second month 
the crossings are spread over a wide arc. 



responding crossing points to result in a loss of output map qual- 
ity due to larger baseline errors. Destriping relies then more on 
the crossing points which are near the ecliptic poles. When the 
observation period is shortened further, only part of the sky is 
covered. Since the crossing point structure is not necessarily 
changed qualitatively, we expect the map quality to worsen more 
slowly as a function of survey duration. Because of the cycloidal 
scanning strategy, the crossing point structure, however, changes 
with a 6-month period, so it will be diff'erent for diff'erent weeks 
or months. 

In Fig. [45] we show the residual / map rms as a function of 
^base for 1 -month, 7-month, and 1-year surveys. The results for 
the Q and U maps look qualitatively the same, except that the 
signal baseline contribution is, of course, much lower. We have 
excluded all pixels with rcond < 0.4 or ^hit ^ 400 from the 
residual maps. 

There are basically two kinds of eff'ects contributing to these 
results. The more trivial eff'ect is that of the lower hit count. The 
average number of hits per pixel in the 7-month survey is 7/12 
of that of the full-year survey, so we expect that alone to in- 
crease the residual map rms by a factor of ^|\2|1 - 1.31. This, 
indeed, accounts for most of the change in the unmodeled 1 // 
contribution to the residual map rms. When the survey segment 
is shortened below 7 months, the number of included pixels falls 
almost in line with the number of samples, so the hit count per 
pixel stays now roughly constant, and there is not much addi- 
tional eff'ect when going from 7 months to 1 month. 

The other eflfect is that of the change in the pattern of cross- 
ing points on solving the baselines. This shows clearly in the 
white noise baseline and signal baseline contributions, where 
baseline errors were important already for the full 1-year survey. 
Here the change from 7 months to 1 month brings also a signifi- 
cant change. A closer inspection of the diff'erent months reveals, 
however, that it is only some of the months that are clearly worse 
than the 7 -month case. 

The first and the second month represent the two extreme 
cases. See Fig. |46l For the first month we are close to the sit- 
uation where all scanning circles cross each other at the same 
pixel. The eff'ect of this is the most striking for the signal base- 
line component, for which it couples to the signal gradients in 
those few pixels where the circles cross and the pattern of hits 
within those pixels. 




1000 



days 



Fig. 47. The binned white noise map rms as a function of sur- 
vey length for ^base = 1 min. Three diff'erent starting points were 
considered: from the beginning of the simulated mission, start- 
ing 30 days later, and starting 61 days later. The black line is for 
the Stokes / map, red for 2, and green for U . 



To study the eff'ect of the length of the survey segment, and 
also its timewise location with respect to the cycloid, we fixed 
the baseline length to ^base = 1 min, and considered survey seg- 
ments of 1, 2, 4, 6, 8, 10, 16, 24, 32, 40, 48, 64, 128, 256, and 366 
days. We also considered the eff'ect of extending the mission to 
488 days (16 months). Moreover, for each survey segment length 
(except the 16 month one) we considered three diff'erent starting 
points for the segment. To separate the eff'ect of the change in 
hit count from the eff'ect of crossing point structure on solving 
baselines, we plot first the rms of the binned white noise map, 
which has only the first eff'ect, in Fig.|47l 

To see the second eff'ect, we consider the ratio of the residual 
map rms to the rms of the binned white noise map from the same 
data. This tells us how well we are doing compared to the white 
noise level. See Fig.|48l We see that for all three starting points, 
the relative map quality consistently improves after 45 days with 
further increase in the survey length up to a full year. The relative 
quality (with respect to binned white noise) of the 16-month map 
is the same as the 12-month map. 

For shorter survey pieces the behavior as a function of sur- 
vey length depends on the starting point. For Q and U , short 
segments, 10 days or less, are clearly much worse than longer 
ones. For the / map, the relative quality in some cases worsens 
with time up to about 40 days. This is related to the changing 
pattern of crossing points in the cycloidal strategy. For some 
periods of time the pattern is more ill-conditioned and adding 
data from such a period to another short period makes things 
worse rather than better. From Fig. [49] we can conclude that one 
such period is somewhere between days 1 and 24 and another 
somewhere between days 70 and 101 of the simulated mission, 
since the residual / map rms is increasing during those periods. 
These coincide with the times when the cycloidal scanning strat- 
egy is producing a clustering of crossing points of nearby scan- 
ning rings at the corners of caustics around ecliptic poles. The 
baselines can be better solved from a more widely distributed 
set of crossing points (Wright 1996). We can see from Fig.l48l 
that the unmodeled 1 // and reference baseline contributions stay 
relatively constant, so the variation indeed comes from baseline 
error. 





1000 



days 

Fig. 48. The ratio of the residual map rms to the binned white 
noise map rms (see Fig. [47]) as a function of survey length for 
^base = 1 min. The black line is for the full residual (excluding 
binned white noise), blue for the white noise baseline map, red 
for residual 1 //, and green for signal baseline map. We show 
also separately the reference baseline contributions {dashed). 
The top panel is for the Stokes parameter /, the bottom panel 
for Q. Results for U are qualitatively the same as for Q. For 2, 
the signal baseline component is plotted multiplied by a factor 
of 100 also {dotted). 



10. Conclusions 

We have described our destriping map-making method 
(Polar/Madam) for CMB surveys. The method has a parame- 
ter, the baseline length ^base, that affects the performance of the 
method. With long baselines, the Madam code is faster and re- 
quires less computer memory. The computer time and memory 
requirements of the code are discussed in Ashdown et al. (.2009 ), 
where it is compared to other codes and methods. 

Here we have done a detailed analysis of the residual errors 
in maps produced with the method. In this paper we concen- 
trated on destriping without a noise prior. For short baselines 
the results can be improved by utilizing prior information on the 
noise spectrum. This will be described in (Keihanen et al. 120091) . 

We have divided the destriping residuals into six compo- 
nents. Three of them, white noise reference baselines, unmod- 
eled 1 // noise, and pixelization noise reference baselines, are 
easy to estimate analytically from the noise power spectrum 



Fig. 49. The top panel of Fig.|48]replotted so that the horizontal 
axis corresponds to the distance of the end of the included data 
from the beginning of the simulated mission, so only for the case 
of the first starting point does it correspond to the length of the 
data used. 



and the signal angular power spectrum. In the map domain the 
baseline components appear as a superposition of thin constant 
stripes of length ^base, whereas the unmodeled 1// noise varies 
along such stripes with a period comparable to ^base- 

The three other components are related to how accurately 
baselines can be solved from crossing points, and depend on 
the scanning strategy. These baseline error components are min- 
imized when there are very many crossing points widely dis- 
tributed. Especially when making partial sky maps from short 
survey segments around a time when ring crossings cluster in 
the same small region of the sky these errors may blow up. Since 
the baseline errors are correlated over long time scales, they pro- 
duce wide bands in the map domain and are important for low 
multipoles. 

The relevance of this analysis is that it guides us in what 
kind of noise residuals to expect in the maps for given detector 
noise spectra, and what baseline length to use in destriping map- 
making. 

For long surveys, with a good distribution of crossing points, 
the dominant residual error components are the white noise base- 
lines and unmodeled 1 // noise. Their combined eff'ect can be 
minimized, when the baseline length is chosen according to 
Eqs. dSSl) or ([6Qb , which put it close to ^base ~ 1/(2/^), where 
fk is the knee frequency of the noise. Because of the other error 
components, one should choose a somewhat longer baseline than 
this. (When a noise prior is used, shorter baselines are better.) 

For a PLANCK-like scanning strategy, where the same circle of 
the sky is observed many times, the diff'erence between baseline 
lengths from the spin period to the repointing period is mainly 
due to nutation, and is small when the nutation is small com- 
pared to the map pixel size. If the knee frequency is comparable 
to, or smaller than, the spin frequency, then the baseline length 
should be chosen from this range, if no noise prior is used. For a 
higher knee frequency the residual errors are larger, and a shorter 
baseline is better. 

Acknowledgements. The work reported in this paper was done by the CTP 
Working Group of the Planck Consortia. Planck is a mission of the European 
Space Agency. We thank K. Gorski, C.R. Lawrence, and J.P. Leahy for useful 
comments. This work was supported by the Academy of Finland grants 205800, 
213984, 214598, 121703, and 121962. We acknowledge the support by the ASI 



30 



Kurki-Suonio et al.: Destriping CMB temperature and polarization maps 



contracts "Planck LFI Activity of Phase E2" and 1/016/07/0 "COFIS". RK is 
supported by the Jenny and Antti Wihuri Foundation. HKS thanks Waldemar 
von Frenckells stiftelse, HKS and TP thank the Magnus Ehrnrooth Foundation, 
and EK and TP thank the Vaisala Foundation for financial support. This work 
was supported by the European Union through the Marie Curie Research and 
Training Network "UniverseNet" (MRTN-CT-2006-035863). We thank CSC 
(Finland) for computational resources. We acknowledge use of the CAMB code 
for the computation of the theoretical CMB angular power spectrum. This work 
has made use of the Planck satellite simulation package (level S), which is as- 
sembled by the Max Planck Institute for Astrophysics Planck Analysis Centre 
(MPAC). Some of the results in this paper have been derived using the HEALPix 
package (Gorski et al. ,2005< . 



References 

Armitage, C, and Wandelt, B.D. 2005, Phys.Rev. D70, 123007. 
Ashdown, M.A.J., Baccigalupi, C, Balbi, A., et al. 2007, A&A, 467, 761. 
Ashdown, M.A.J., Baccigalupi, C, Balbi, A., et al. 2007, A&A, 471, 361. 
Ashdown, M.A.J., Baccigalupi, C, Bartlett, J.G., et al. 2009, A&A, 493, 753 
Borrill, J. 1999, in Proceedings of the 5th European SGI/Gray MPP Workshop, 

Bologna, Italy, astro-ph/99 11389 
Burigana C, Seiff'ert M., Mandolesi N., Bersanelli M. 1997, Int. Rep. 

TeSRE/CNR, 186/1997, March 
Burigana, C, Malaspina, M., Mandolesi, N., et al. 1997, Int. Rep. TeSRE/CNR, 

198/1997, November, | astro-ph/9906360| 
Burigana, C. & Saez, D. 2003, A&A, 409, 423. 

Challinor, A., Fosalba, R, Mortlock, D., et al. 2000, Phys.Rev.D, 62, 123002 
Couchot, R, Delabrouille, J., Kaplan, J., & Revenu, B., 1999, A&A Suppl. Ser 
135, 579 

de Gasperis, G., Balbi, A., Cabella, R, NatoH, R, & Vittorio, N. 2005, A&A, 

436, 1159. 
Delabrouille, J. 1998,A&ASS, 127, 555 

Dore, O., Teyssier, R., Bouchet, F.R., Vibert, D., & Prunet, S. 2001, A&A, 374, 
358. 

Dupac, X., & Tauber, J. 2005, A&A, 430, 363 
Eftstathiou, G., 2005, MNRAS, 356, 1549 
Eftstathiou, G., 2007, MNRAS, 380, 1621 

Ganga, K. 1994, Ph.D. Thesis, http://spideripac.caltech edu/staff/kmg/ 
Goldberg, J.N., Macfarlane, A.J., Newman, E.T., Rohrlich, F., & Sudarshan, 

E.C.G. 1967, J. Math. Phys. (N.Y.) 8, 2155 
Gorski, K.M., Wandelt, B.D., Hivon, E., Hansen, F.K., & Banday, A.J. 2005, 

ApJ, 622, 759 

Harrison, D.L., van Leeuwen, F., & Ashdown, M.A.J. 2008, in preparation 
Keihanen, E., Kurki-Suonio, H., Poutanen, T., Maino, D., & Burigana, C. 2004, 
A&A 428, 287 

Keihanen, E., Kurki-Suonio, H., & Poutanen, T. 2005, MNRAS, 360, 390 
Keihanen, E., et al. 2009, in preparation 
Keskitalo, R., et al. 2009, in preparation 
Leahy, J.P., et al. 2009, in preparation 

Maino, D., Burigana, C, Maltoni, M., et al. 1999, A&ASS, 140, 383 
Maino, D., Burigana, C, Gorski, K.M., Mandolesi, N., & BersaneUi, M. 2002, 
A&A, 387, 356 

NatoH, R, de Gasperis, G., Gheller, C, & Vittorio, N. 2001, A&A, 372, 346. 

Newman, E.T. & Penrose, R. 1966, J. Math. Phys. (N.Y.) 7, 863 

Planck Collaboration 2005, The Scientific Programme of Planck ("Planck 

Bluebook'), ESA-SCI, 2005, 1 ( astro-ph/0604069 I 
Poutanen, T., Maino, D., Kurki-Suonio, H., Keihanen, E., & Hivon, E. 2004, 

MNRAS, 353, 43 
Poutanen, T., de Gasperis, G., Hivon, E., et al. 2006, A&A, 449, 1311 
Press, W.H., Teukolsky, S.A., Wetterling, W.T., & Flannery, B.R 1992, 

Numerical Recipes, 2nd ed. (Cambridge University Press, Cambridge) 
Reinecke, M., Dolag, K., Hell, R., Bartelmann, M., & Ensslin, T. 2006, A&A, 

445, 373. 

Revenu, B., Kim., A., Ansari, et al. 2000, A&ASS, 142, 499 
Sbarra, C, Carretti, E., Cortiglioni, S. et al. 2003, A&A, 401, 1215. 
Seiff'ert M., Mennella A., Burigana C, et al. 2002, A&A, 391, 1185 
Stompor, R. & White, M. 2004, A&A 419, 783 

Tuovinen, J. 2003 Planck Newsletter Number 4, p. 7 (www.rssd.esa. 

int/SA/PLANCK/docs/Newsletters/PlanckNewsletter4.pdf) 
van Leeuwen, F., Challinor, A.D., Mortlock, D.J., et al. 2002, MNRAS, 331, 975 
Varshalovich, D.A., Moskalev, A.N., & Khersonskii, V.K. 1988, Quantum 

Theory of Angular Momentum (World Scientific, Singapore) 
Wandelt, B.D. & Gorski, K.M. 2001, Phys. Rev. D63, 123002 
Wright, E.L. 1996, Report UCLA-ASTRO-ELW-96-03, November, 

astro-ph/96 12006 
Yvon, D. & Mayet, F. 2005, A&A, 436, 729 
Zaldarriaga, M. & Seljak, U. 1997, Phys. Rev. D55, 1830 



